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“ AESTRACT 


The two principal existing methods of calculating axially- 
symmetric compressible flow in — are: (1) a simplified 
one-dimensional analysis, and (2) numerical methods using the com- 
plete or linearized flow equations. The first is not satisfactory for 
multi-stage turbines with appreciable wail divergence; the second is 
very tedious and time consuming. The purpose of this investigation 
is to extend the approximate methods, successfully used in calcu- 
lating incompressible flow in compressors with constant blade height, 
to the analysis of compressible flow in turbomachines with variable 
blade height. Assuming that the blades can be completely defined by 
the exit flow angle, and neglecting the influence cf downstream blades, 
the analysis is made considering the flow between successive blade 
rows only. With these restrictions, subsonic and isentropic super- 
sonic flow patterns can be determined for arbitrary boundary shapes 
as long as separation does not occur. Average losses can be ac- 
counted for by the use of a polytropic law, and the effect of radial 
variations in stagnation temperature can be included without diffi- 
culty. Examples illustrating the flexibility and practical value of the 
iteration method, and the rapid convergence of successive solutions 


arë given. 





are 
I. INTRODUCTION 

In the modern aircraft gas turbine, the flow is necessarily 
in the transonic range for two reasons: (1) to keep the weight and 
size of the machine at a minimum, and (2) to fit tre mass flow and 
pressure ratio of the compressor. Earlier turbojet engines had single 
stage transonic turbines, but present engines employ multistage tur- 
bines with rapid channel divergence in order to use effectively the 
increased pressure ratios and keep the axial velocities compatible 
with the blade speed. In these types of machines, the flow at the 
mean blade height is mainly axial, and the aspect ratio of the blades 

P a usually above unity for weight saving considerations. As long as 
the channel wall divergence is not too great so that separation can 
be avoided, there is no reason for low efficiency, particulariy in 
view of the favorable pressure gradient through the turbine. 

Two distinct design problems occur in practice which Marble 
(Reference 15) has classified as the direct and inverse problems of 
turbomachine design in analogy with the corresponding classical 
problems in the theory of finite wings. The direct, or '"oíf-design' 
problem occurs when the blade speed, blade shape, boundary coníigu- 
ration, and fluid state ahead of all blades is given, and the threc- 
dimensional velocity field, blade loading, and distribution of energy 
are to be determined. The inverse problem occurs in the original 
design of turbomachinery, when the blade loading, blade speed, 
boundary configuration, and fluid state ahead of all blades is pre- 
scribed, and the three-dimensional velocity field, blade shape, and 
distribution of energy in the field are to be determined. 


The exact analysis of viscous, compressible flow through an 
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axial turbomachine will probably never be possible because of the 
non-linear character of the cquations governing the flow. However, 
the design engineer has a great need for an analysis suitable for 
engineering purposes, and may have no use for a lengthy, compli- 
cated and costly procedure even if it were available. He needs an 
extremely flexible method of analysis which, while being as simpli- 
fied as practicable, still produces inforraation accurate enough for 
effective design work. 

the equations governing three-dimensional, cornpressidle 
channel flow are greatly simplified if terms involving viscosity are 
neglected. In view of the fact that viscous effects actually are very 
slight except near the boundaries, and as long as separation is avoided, 
the assumption of zero viscosity seems justified. Viscosity efiects 
could be accounted for as an added correction to the basic solution. 
However, with the present day turbines which operate in the high 
subsonic range, it is obvious that compressibility effects can not 
be neglected. 

One of the principal methods available to the design engineer 
for the analysis of axially-syrmmetric, compressible flow in turbo- 
machines is a simplified one-dimensional analysis. This method is 
used extensively in preliminary design considerations and is simple 
and straightforward. Since only the exit flow angles from the blade 
rows at the mean blade height are considered, the radial distributions 
of velocity and energy are not determined. Although this analyeis 
gives reasonably accurate results for a single stage turbine witha 


fairly high hub ratio and small channel divergence, the one-dimensional 


method is inadequate for large channel divergence where the radial 


23% 
distribution of velocities becomes significant. It is particularly in- 
adequate for the design of multistage turbines with channel wall di- 
vergence since it is difficult to get the prover balance of worl on the 
individual stages. This work balance is strongly affected by the radial 
distribution of the axial and radial velocities. Since the one-dimen- 
sional analysis is done at the mean blade height, it gives no informa- 
tion as to the radial distribution of velocities. 

The second principal method of analysis of axially-syimmetric 
compressible flow in turbomachines is by numerical methods using the 
complete or linearized flow equations to solve for the complete velocity 
and energy field through the machine. Since the determination of the 
entire flow field involves six unknowns at each point in the flow (the 
three components of velocity, the pressure, the temperature, and the 
density), the analysis using these numerical methods is necessarily 
lengthy and extremely complicated and cannot be considered suitable 
from the point of view of the design engineer. Monroe (Reference 7) 
formulated the idealized problem in terms of a stream function for the 
velocities in the meridional plane, and solved the resulting non-linear 
differential equation by a simultaneous application of the relaxation 
technique of Southwell and an iteration process. Wu (Reference 9) has 
utilized various relaxation methods, matrix systems, and finite differ- 
ence schemes to solve the complete hydrodynamical equations for the 
entire flow field in many types of problems. Vavra (Reference 8), 
Reissner (Reference 10), Goldstein (Reference 11), Wislicenus (Ref- 
erence 13), and many others have employed various techniques to solve 
the complete, or nearly complete, equations by some numerical method. 


Although these numerical methods are not suitable for design purposes, 
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they are quite general and can be used for centrifugal, mixed ilow or 
axial flow machines with blades of any aspect ratio. 

In an effort to simplify the problem of solving for the complete 
flow field, a linearized analysis has been provided by Marble (Refer- 
ences 5 and 15) and Fabri (Reference 6). Marble (Reference 15) pro- 
posed a linearized analysis of the problem, satisfactory for high aspect 
ratio blades, to allow an approximate treatment of the general blade row 
with prescribed wall geometry. He assumed that the trailing vorticity 
could be considered to be transported downstream by the mean axial 
flow and not by the perturbation velocities, that is, disturbances in the 
radial and axial velocities are small in comparison with the mean axial 
velocity, and that the blade row is made up of an infinite number of 
blades of finite chord. The analysis for the inverse problem was carried 
out in detail and was found to allow a general and reasonably simple so- 
lution. In Reference 5, Marble extended the analysis to the direct prob- 
lem of off-design operation, and to the study of mutual interference of 
neishboring blade rows in a multistage axial turbomachine. He provided 
examples of axial and conical flow, and stated that the simple linearized 
solution was sufficiently accurate if the vorticity effects were not large. 
A second order linearization was given to handle problems with greater 
vorticity effects. However, the computation required to find the three- 
dimensional flow in any particular case is rather lengthy using the lin- 
earized solutions, and involves some type of numerical integration. A 
simple exponential approximation was then introduced to simplify the 
computation, and the results using this approximation compared favor- 
ably with the detailed results of the linearized solution. 

The analysis of flow in turbomachinery where the channel bound- 


aries are purely cylindrical, or vary only slightly from a cylindrical 


A 
shape, has becn extensively investigated by Rannie (References 1 and 14), 
Traupel (References 2 and 3), Sinnette (References 4 and 18), and others. 


In these analyses, the radial momentun: equation used in the treatment is 


of the form, 


ven, We 
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o Ər (1), 
where the radial force and other moinentum terms have been ncglected. 
The distribution of axial velocity is approximated by assuming axial sym- 
metry, and since the centrifugal force within the rotating fluid body is bal- 
anced only by the radial pressure gradient, the resulting u calculated 
in this manner is in reality that which must exist far downstream of the 
blades where radial velocities and accelerations have vanished. No in- 
formation is provided by these analyses on how rapidly the change in ve- 
locity pattern takes place as the fluid passes through the blade row. 

The purpose of this thesis is to present an approximate method 
for the analysis of compressible flow in modern axial turbomachines which 
gives sufficiently accurate information for design purposes and yet retains 
the simplicity necessary to be of practical value to the design engineer. 
This is the logical step forward from the one-dimensional method to the 
analysis of compressible flow in axial turbomachines. Methods which 
have been successfully employed in the incompressible case (Reference 1), 
have been modified and extended to the analysis of compressible flow in 
arbitrarily shaped azially-symmetric channels. It is obvious that little 
simplification over the so-called "numerical" methods can be obtained if 
the complete flow and energy fields are to be determined in the analysis. 
In this thesis, the flow is analyzed at the station immediately downstream 
of the blade row by a consideration of the upstream conditions, Channel con- 


figuration, and blade shapes. The blade trailing edge angle only is consid- 


ered, and the complete flow field through the blade row is not determined. 
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The method of analysis presented is extremely flexible and kas 
been used to solve the direct and inverse problems, subsonic or Super- 
sonic, for arbitrary axial boundary configurations. In the solution of 
the direct problem, the choice of blade speed, blade spouting angle, 
boundary configuration, and upstream conditions car be arbitrary, and 
the threc-dimensional velocity field, blade loading, and distribution of 
energy in the field are determined. In the inverse problem, for a given 
blade loading, blade speed, boundary configuration, and fluid state ahead 
of all blades, the blade spouting angie, velocity field and energy distri- 
bution can be determined. By a simple change of sign in the equations 
used, the method can be applied equally well to an axial compressor as 
to an axial turbine. Aithough this method is not as general as those used 
by Wu (Reference 9), for example, it is more ilexible and practical for 
blading oí aspect ratio greater than unity, and is sufficiently accurate 
in view of the assumptions used in the treatment. 

The analysis is carricd out on the basis of steady, axially- 
symmetric flow with infinitely many blades. The fluid is assumed to 
be compressible, but all other real fluid effects are neglected except 
as they may be approximated by a "'polytropic efficiency'' and the use 
of the polytropic exponent ''n''. As a result of this assumption, the pos- 
sibility of separation of flow at the channel walls or along the bladcs is 
excluded. Boundary layer effects are also not considered in the anal- 
ysis. In the supersonic solution, the occurrence of shock waves is 
excluded. Radial or nearly radial blades are presupposed so that the 
radial force can be neglected, and the radial momentum equation used 


in the treatment has the form, 


4 0D _ Es PU O0 (2), 
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This assumption Rei justified in view of the investigations oí Karlsson 
(Reference 19), who, by a two-dimensional analysis of incompressible 
flow with an infinite number of blades, showed that the cffect of the radial 
force resulting from normal blade twist is of negligible magnitude. 

The effect of downstream blades on conditions upstream is nc- 
glected, and the blade shape is specified by the trailing edge angle only. 
Simply stated, a blade row is assumed whose only effect on the flow is 
to turn it through a specified angle distribution which is a function of 
radius across the T An assumption similar to this has been used 
successfully in the analysis of flow through a compressor (References 
1 and 14), and should be even more applicable to turbines where the blade 
solidity is high. The cquations are particularly suited to the analysis 


of flows «with a prescribed radia! total enthalpy distribution. 
R 


v 


Subject to the above assumptions, the resulting "exact" equations 
are solved by an approximate method using a simple iteration procedure. 
An integral for the axial velocity is derived in terms of known upstream 
conditions and the downstream pressure. The pressure distribution is 
found by integrating the radial momentum equation and applying the 
energy equation at a boundary. The mass flow equation is used to solve 
for the constant obtained from the axial velocity integral. The itcration 
is started by assuming an axial velocity distribution and continued until 
the solution is reached, i.e., until the velocity distribution produced by 


an iteration is the same as that assumed for that particular itcration. 





Be 


Iii. FUNDAMENTAL EQUATIONS 


A. Development of Basic Differential Equations 


The differential equations are developed assun.inz steady, 
compressible flow of a fluid through an axially symmetric cylindrical 
channel of arbitrary shape. An infinite number of blades is assumed 
so that the velocity components are considered independent of angle 
about the axis of the machine. The fluid is assumed perfect with no 
losses of any sort except as these losses can be approximated by a 
"polyíropic efficiency" and the use of the polytropic exponent "n", 

All boundary layer and viscosity effects are also neglected. 

In order to generalize the equations for an arbitrary (but 
reasonable) channel shape, it is necessary to consider a special type 
of blading. These blades, infinite in number, have no other effect on 
the flow than to turn it through a prescribed angle distribution. Thus 
the effect of downstream blades on upstream conditions is neglected. 
Attention is directed to the flow patterns immediately before and after 
the blade rows for either the stator or the rotor. 

In Figure |, a stream annulus of small radial extent is shown, 
and absolute velocity components are considered at stations 1, 2, and 
3. Generally, there will be a shift of the stream surface in passing 
through either stator or rotor blading, and three components of 
velocity willbe present at each station. 

With the above assumptions, the continuity equation is written 


in differential form as: 


PMT dy = pru dz o eun a), 
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BLADING 






STATION | STATION 2 STATION 3 


FIG. 1 ——- CHANNEL AND STREAM ANNULUS 

Tor a row of stator blades, the energy equation must hold along 
a stream surface, and for a row of rotor blades, the energy equation 
must hold along a relative stream surface. Hence, in terms of 


absolute velocities at stations 1, 2, and 3, we have: 


<a 1B f UFV + *) = x zu (uti E (4) 
x Fe / "e 2 2 i = y 79% / 7 e 


. tadial force, the radial momentum equation must 
hold at all stations. In absolute velocities and omitting subscripis, 


this equation is: 


star . eee. ae (2) 
e a7 p» ^s Waz : 


If the flow is considered to be isentropic along 2 relative stream- 


surface, the following relations must hold along this surface: 





iQ = 


E D 
Um . PER TEN t 6 
ce (BY: GF) (6) 
I lt: 
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If a polytropic flow is considered, 7 will be replaced by n in 
Equations 6 and 7. The above equations are fundamental in the solution 


of the direct aui inverse problems. 


B. Transformation to Non-dimensional Form 


For ease in handling the equations, and to avoid the use of large 
numbers in calculation, the derivation is completed in non-dimensional 
form. The notation used is given in Appendix A. Axial and radial cc- 
ordinates are made non-dimensional by referring them to a standard 
dimension, usually the outer radius of the annulus at station i. Velo- 
aties are made dimensionless by referring them to a standard tip speed. 


The following notation is used: 

















R - outer radius ci the annulus at station |, 
o 
A Ro - reference blade speed. : 
a e A . e. 
d- w fe - dimensionless radial cornponent of velocity. 
v i ; 
À= OR - dimensionless tanpential component of velocity. 
o 
$- Ey - dimensionless axial component of velocity. 
eo 
Te. - dimensionless pressure. 
Po 
EDS : : 
Q = ° e - dimensionless density. 
X 
R ' 
T= 22d - dimensionless temperature. 
wer; 
a = 7 - gas equation. 





wil 





m = We - dirnensionless mass flow, 
2719 72 Re 
é 2 X - dimensionless radial distance. 
o 
M 5 ~ dimensionless axial distance. 
With this notation, Equations (3), (4), (5), (2), (6), and (7) become 
respectively: 


u" eic A A T 


C. Equations for the Direct Problem for an Arbitrary Channel Shape 


The basic integral for the axial velocity will be derived for the 
rotor (stations 2 to 3) for.the isentropic case and simplified to the 
integral for the stator. The corsplete derivation is given in Appendix 
B for the isentropic direct problem and in Appendix C for the polytropic 
direct problem. | 


Differentiating the left side cf Zquation (10) in respect to the 


radius 3 and the right side in respect to £ gives: 
3 





ar_ at Ber. A), 1 & [gh 2, 42), AGA) 
N Far: "hag u), 


Using Equations (8) and (13), differentiation of the first term of the right 


side of the above equations results in: 


7 d 2) E dad, 7 2d 7 ds, (15), 


7| d£N9/  Qd£, +o Hat E 


Substituting Equations (11) and (15) into (14), and using the differential 


d £e 





form of the mass flow relation (i2quation (8) ) for the ratio 


| | d $ 
the following differential equation is obtained: 
AG 
A, eM. «d 2 Pa Z £5] 7 oler al$2A2) 
£ z ag, Pa) Peña | 71 af 9£, 
(16) 


P. | 


For a cascade with an infinite number of blades, the exit flow 
is parallel to the tangent of the camber line at the trailing edge, and 
its direction is independent of the direction of flow approaching the 
cascade. For most practical applications, the direction of flow 
leaving the cascade is very nearly constant and sr of the 
inlet direction throughout the design operating range. With this 
assumption for the rotor, it is therefore assumed that the rotor 
spouting angle Ps is a known function of radius at station 3. Using 
the sign convention shown in Figure 8, ii is therefore possible to 
describe the relation between Az and Pz along the radius at 


station 3 with a function G ( £,) such that: 


= GE )B- £, an, 
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Eliminating A from the differential Fauation (16) by use of 
Equation (17) and integrating, the following integral for the rotor axial 


velocity for the direct problem is obtained: 
fia [6 
E (I+ > 
O E 2E dii. d (fade) 
3 £st NET p £ 2-, de d£ 
f; 











- 250) las A), yaar (erie iar 
(18), 

It is to be noted that no linearization has been performed in 
arriving at the above integral. The expression for $3 is quite cornpli- 
cated, but its numerical evaluation by 2 simple iteration process to be 
described later is straightforward. Once a set of streamlincsa have 
been assumed, all factors involving the upstreain conditions becorne 
known functicns of 5 3 « Inthe solution of the direct problem, it 
is assumed that C (83) is a given function of 7 . Thus in the 


integrand, P3 and on are the unknowns which must be determined 
3 





in the iteration process in order to evaluate the integral. 
An alternate form for the square bracket in the integrand can be 


obtained by using the energy and radial momentum equations at station 2: 


y BIS +S = ¿qe +&,)- i s) zg fA Jal 1) (19), 


In the case that P,, = l (isentropic expansion to station 2 from stagna- 
tion conditions), the bracketed terra in the integrand may be further 


simplified to: 


2 (7) 6) e 





Jas 


The constant of integration Cy is determined by integrating 
Equation (8). The integrated form of this relation equates the total 


flow rate upstream to that at station 3: 


PR 
m 7 Q, 22) A (21). 


Using Equation (13), the integral form of the radial momentum 


equation becomes: 


Tz 
BF (E) a (2 ATA a (22). 


The constant of integration C3 is obtained by applying the energy 





relation (Equation (10) ) at a boundary. 


Having assumed a set of streamsurfaces for a particular 





de; i 
can be approximated 
d $s - 


from the values of the slope and rate of change of slope of these stream- 


boundary configuration, the value of 


surfaces at station 3. At each point on a streamline at station 3, the 


following relation must hold: 











$ E 3 (23), 
where GI J =), is the slope of the streamline in question, Thus 
the value of je becomes: 

2 
ip i = (EE) (24), 
In order to approximate the value of qa at points across the 


channel, it is therefore necessary to estimate the slope and curvature 


of the assumed strearnlines at each point. 
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For the solution of the direct problem for the rotor with an 
arbitrary channel shape, Equations (10), (13), (17), (18), (21), (22), 
and (24) are applicable, 

The derivation of the integral for the axial velocity in the case 
of the stator is carried out in a similar fashion using the form of the 
energy relation given as Equation (9). Using the aign convention of 


Figure 8, Equation (17) for the stator (stations 1 to 2) becomes: 
A2=F(,) pe (25). 


The integral ior Po obviously takes the same form as 


Equation (13), and is given here only for completeness: 

thds "126 TE 

h = Jdem f^) , PY £e | v atte 
e me 14 PELI AZ [7 2, 
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For the solution of the direct problem for the stator with an 
arbitrary channel shape, Equations (9), (12), (21), (22), (24), (25), 


and (26) (with appropriate subscripts) are applicable. 


D. Equations for the Inverse Problem for an Arbitrary Channei Shape 


In the inverse problem for the rotor, the same assumptions and 
basic differential equations are used as in the direct problera with the 
exception that now the blade relation, Equation (17), is replaced by 
the following given function which describes the work ontout of the 


rotor as 53 function of the radius at station 3: 


FG) = ace ee (27), 


t= 


The function H /( N) gives the dimensionless rate of change 
of angular momentum per unii mess flow at the particular radius, and 
represents the energy taken out of the stream. By multiplying this 
function boy suitable constants it can readily be converted to D. T.U. 
per pound of mass flow, 

The complete derivation of the axiai velocity integral is given 
in Appendix D for the isentropic case and in Appendix E for the poly- 
tropic process. For the isentropic inverse problem, the derivation 
remains the same up to and including Equation (16), which may be 
written in a siightly modified form as: 


3 BZ Y et 
2, tag led g) 2 az fs A;) = (2 hfe a- age.) 


Bs 
2-4 7 
= RL (E + h LZ 
287 2) PE (28). 
Using Equation (27) to eliminate A, , tre integral is 





5 
e 7 
p = f £s iz ales _ xA ug (ai) (H - SEE Rea] d d 65 


Ai $a 7 dZ, £ 
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Ps a 7 (29). 


Again the numerical evaluation of this integral by the iteration 


— process to be described later is straightforward once a set of stream- 
d 83 


lines has been assumed, The unknown P, and a£ appear in 
3 





the integrand as in the direct problem. Because vi the nature of the 


ya 


differential equation, however, the axial velocity C3 algo appears in 


the integrand., 


zs 


Since the mass flow and radial momenturi intezrals remain 
unchanged, Equations (10), (13), (21), (22), (24), (27), and (29) are 
applicable to the solution of the isentropic inverse problem. 

If the inverse problem could be conceivably ER. to a stator 
it would consist of prescribing values of A 2 as a function cf A , 
and solving for the resulting velocity field, distribution of energy and 
blade shape (as described by the blade trailing edge angle). In this 


case, the integral for the axial velocity simplifies to: 


2 
£ E 
b p n E ates 
fi (30) 
where Aaló,) is a prescribed function of the radius at station 2. 
Equations (9), (12), (21), (22), (24), and (30) (with appropriate 


subscripts) would be applicable. 


E. Changes in Equations for the Compressor Solution 


In the development of the equations in the preceding sections, 
the turbine sign convention given in Figure 8 has been adopted. If 
the equations are to be used for compressor design, Equation (17) for 


the direct rotor problem becomes, in accordance with Figure 9: 


= (31), 


and the energy relation, Equation (10), is changed to: 


A L(+ i+ BF -afA) zZ 1g x -2£2,) (32). 


7-1 


~~ We 


Using these equations, the derivation is carried out in exactly the 
same manner as set forth in the preceding sections for the turbine. 
For the inverse problern for the compressor, the work out- 


put function (Equation (27) ) becomes: 


H( )- fe Je 425 


Using Equations (32) and (33) in piace of Equations (10) and (27), the 


(33). 


relations for the inverse compressor problem are derived in the same 


manner as for the turbine. 
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Hi, PROPOSED ITERATION PROCEDURE 
A, Procedure for the Direct Problem 


The general scheme of iteration is exactly the same for the 
direct and inverse problems, subsonic or supersonic, isentropic or 
polytropic. In Appendix F, one complete iteration is given for an 
isentropic supersonic solution (Example Vb) and a polytropic subsonic 
solution (Example VIII b) for the direct rotor problem. The stator solu- 
tion is obtained by the same procedure using the appropriate axial 
velocity integral and energy equation. 

In order to understand the technique of iteration for the direct 
problem described below, the reader is urged to refer frequently to 
axamples Vb aad VIH b of Appendix Y. for an arbitrary channel 
shape, seven equations are used in the solution. For the rotor, these 
relations are Equations (10), (13), (17), (18), (21), (22), and (24). In 
the examples worked out in Appendix F, the functions Y (£, ) and G (5) 
have been assumed to be constant and independent of radius in order 
to simplify the axial velocity intecral and the computation. This 
assumption is by no means necessary. For arbitrary given functions, 
F ( ee) and G ( =) , the exponential integrals of Equation (18) 
should be evaluated at various stations across the channel by plotting 
the integrand and determining the value of the integral graphically. 

In the solution of the direct problem, it is assumed that the 
boundary configuration, blade speed, blade trailing edge angle, and 
upstream conditions are prescribed. It is further assumea that the 
mass flow rate is fixed and is compatible with the given upstream 


conditions. 


¿Ae 


e. 
Once the problem has been set up, the first step is tc assume 


a set of streamsurfaces as shown in Figures 3, 4, and 5. For most 
applications, these streamsurfaces should be chosen on the basis of 
equal annular areas in a plane perpendicular to the axis of the machine. 
The accuracy of the solution can be improved by increasing the nuraber 
of stations across the channei. However, itis believed that for most 
purposes, six or seven steps will give sufficient accuracy for design 
work. 

Having assumed a reasonable set of streamlines, the second 
step is to estimate the axial velocity distribution at the downstream 
station. Taking average or middle-channel values for the prescribed 
upstream conditions, the mass flow relation and the energy equation 
can be cornbined into one equation in terms oi $3: Using average 


ya 


values in Equations (10) and (21), the relation is obtained as follows: 


T= Toe- Pm rH EXHEX) gn 


d 


m- gaa EE) (35). 


Substituting Equation (34) into Equation (35) gives: 
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Using Equation (17) to eliminate Az and Equation (23) to 
eliminate O, , Equation (36) becomes a relation involving only 
one unknown, $; . The solution must be found by trial and error. 


In general, there will be two values of $5 which satiafy this 
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FIG. 2 — PLOT OF MASS FLOW VS. AVERAGE AXIAL 
VELOCITY FOR CHANNEL FLOW 

If no value of $, satisfies Equation (36), the prescribed 
mass flow cannot be forced through the channel because of choking, 
and the boundary configuration or upstream conditions must be altered 
accordingly. 

Waving found the average value of the axial velocity at station 3, 
an estimate of the distribution of this velocity across the channel can 
be made considering the total temperature gradient, the blade shape, 7 
and the pressure distribution expected at the downstream station. 

Using the estimated distribution of $3 , Equation (22) is 
iterated across the channel to determine the pressure gradient. The 
tangential velocity is found from Equation (17), the radial velocity from 
mquation (23), and the value of ge 


3 
: a Os AO; 
to evaluate E ; az. , and dE, 





from ¿quation (24). In order 








at various points across the 





Be 


channel, itis necessary to estimate as accurately as possible the slope 
and curvature of each streamline at the station in question. Xnowing 
the complete channel configuration, slopes and curvatures at points 
across the channel can be easily estimated using a simple linear distri- 
bution scheme, or a large scale drawing can be constructed and the 


values of slope and curvature taken graphically from this. The value 


d. $5 
d. 75 


stations along the axis. However, it is to be noted that even a rela- 





of in Equation (24) can be taken as the average value between 
tively large error in the estimation of the above quantities wil) have little 
effect upon the pressure gradient obtained because of the nature oí the 
integrated form of Equation (22). The actual pressure distribution 
across the channel is found by applying the energy equation at a boundary 


to determine the constant C 2* 

d. Os 
d E 
integral can now be iterated across the channel for values of $3 in 





With these values of P} and » the axial velocity 


terms of the constant C. If the curvature of the channel is not large, 








the value of E is very small compared to the other factors in 
3 
the integrand (2G, for example), and hence a relatively large error in 
estimating ge has a small effect on the resulting value of E 
3 


The constant C, is obtained by iterating the mass flow relation 
(Equation (21) ) across the channel, and a new axial velocity distribution 
ig obtained. 

The above described process is the procedure fcr one comnlete 
iteration of the given equations. Normally in the subsonic case, the 
distribution of $3 thus obtained is used in the next iteration, and the 


process is repeated until the axial velocity distribution produced by an 


23: 


. 
iteration is the same as that assumed for that particular iteration. 
This is the solution. Values of L; and Q4 are then found by 


use of Equation (13). 


b. Procedure for the Inverse Problem 


ee 


Generally, the same procedure of iteration, as described in 
section 4 above for the direct problem, is used in the solution of the 
inverse problem. Equations (10), (13), (21), (22), (24), (27), and (29) 
are applicable to the isentropic rotor. Relations used in the poly- 
tropic case are derived in Appendix E. Examples IX, X, anu si of 
Appendix F are typical solutions to the inverse problem. One complete 
iteration is given in Example IX b to illustrate the method used. 

It is noted that the integrand of Equation (29) contains a term 
involving A . Since in any practical solution the axial velocity 
at any point across the channel will not be zero, this term does not 
cause any real difficulty in tne iteration procedure. ilowever, if the 
original estimate of the axial velocity profile at the downstream Station 
differs greatly from the actual solution, it has been found that succes- 
sive iterations vill "oscillate'" about the real solution, but will 
eventually converge as in the direct problem. 

it is obvious that the work output function I: ( £ D , given as 
Equation (27), must be prescribed reasonably. For any given channel 
shape, blade speed, upstream conditions and mass flow, there is a 
limited range of blade angles that can be used before the smooth flow 
over the blades is destroyed by separation effects. The limitation thus 


imposed by the maximum turning angle through the blade row also 


= 


Dim 


limits not only the amount of work that can be extracted, but also the 


distribution of this work across the channel. 


C. Converpence of Subsonic and Supersonic Solutions 


Figure 2 indicates that, for any given setup, there is normally 
a supersonic as well as a subsonic solution for the direct and inverse 
problems. This is illustrated in Appendix Y by Examples IVa and Va 
for the stator and Examples IV b and V b for the rotor. In these 
exainples, for a given boundary configuration, blade spouting angle, 
and set of upstream conditions, both the subsonic and supersonic 
solutions have been obtained. : 

Experience in working various examples has shown that if the 
flow is in the low subsonic range, point A of Figure 2, for example, 
the iteration process converges very rapidly, and only two or three 
iterations are necessary to obtain accuracy of four significant figures 
in the solution. However as the flow becomes highly subsonic as at 
point B, the convergence becomes less rapid, and more iterations are 
necessary to obtain the solution. It is interesting to note that, in the 
subsonic range, the process will converge from any reasonable 
assumed distribution of $3. it is even possible to obtain a solution in 
which the flow is supersonic in one part of tae channel and subsonic in 
the remainder (Examples la and Ula). 

It has been found, however, that if the flow is completely 
supersonic (point E of Figure 2, for example), the iteration process 
diverges rapidly away from the solution, Thus if an axial velocity 
distribution is assumed whose average value is represented by point D, 


successive iterations will result in decreasing values of $4 and 
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eventually the ne solution at point A will be reached. If an 
axial velocity distribution corresponding to point F is assumed, 
successive iterations will rapidly increase $3 toward infinity. It 

is the realization of this phenomenon which makes it possible to 
obtain the isentropic supersonic solution in a given problem. Having 
made two successive iterations at each of the points D and F, and 
having noted the relative change in the constant C; in these iterations, 
a simple interpolation will give the values of $4 corresponding to the 
solution at point E once the solution has been ''straddled". 

It is important to understand the basic reason why the proposed 
iteration process is successful. For either the direct or inverse 
problems, six unknowns must be solved for at the downstream station. 
These are the three components of velocity ( 9, , A, and $3) the 


pressure Pa, the density Q}, and the temperature 2, . Using 


the isentropic (or polytropic) relation along streamsurfaces, two of 


% 


these (93 and ur ) cen be eliminated. The remaining four unknowns 
are made functions of the axial velocity $3 and the boundary and 
streamline configuration. The radial velocity is eliminated in terms 
of ¢, by Equation (23), the tangential velocity A, by Equation (17) 
or (27), and the pressure P4 by the integral, Equation (22). Thus, 
effectively, only one unknown remains, the axial velocity 63 5 and for 
any assumed distribution of ¿$ 3° all other unknowns are determined. 
An iteration process must be used to determine $4 because cí the 


coraplexity of the equations involved. 
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D. Correction for the Displacement of Streamlines 

rlaving assumed a set of streamsuriaces for any particular 
example, and having found the corresponding solution by the proposed 
iteration procedure, it is possible to make a correction to determine 
more accurately the actual position of the streamsurfaces and thus 
refine the solution. The mass flow functions m (£) and m (%) 
are plotted against radius for conditions at stations 2 and 3 as shown 
in Figure 6. Using this diagram and fixing the radial coordinates at 
station 3, it is possible to estimate new corresponding upstream 


coordinates. With these corrected upstream radial coordinates, and 


the upstream conditions corresponding to them, the problem is re- 


worked and a new solution found. itis to be emphasized that this is 


a small correction and that the new solution will not vary more than 
one or two percent from the solution based on the oripinal streara- 


surfaces. An illustration of this method is given in Example V b of 


Appendix F. 








ey. 
IV. ADDITIONAL CORRECTIONS AND SIMPLIFICA TIONS 


A. Corrections to tne Basic Solution 


If it is desired to solve for the flow pattern and distribution of 
energy in the field at any axial station in the blade row itself, this can 
be done approximately by multiplying the axial velocity by a function 
which accounts for the blade thickness and thus reduces the effective 
area of the channel. Assuming a blade thickness function t ($) and 
a blade spacing s, the increased axial velocity at any point within the 


blade row would be: 
(4), = ser) Y 


where s and t (€) are assumed to be known functions at any racial 
position across the channel. The assumptions used in the direct 
problem, including that of axial symmetry, remain unchanged, except 
that the blade spouting angle now becomes the mean cainber line at 
the point in question. The fundamental equations and iteration "idm 
cedure remain unchanged except for the use of the above indicated 
expression for the axial velocity. 

One of the basic assumptions used in the derivation of the 
fundamental equations given in Part II, is that the effects of the vorti- 
city of the flow downstream on upstream conditions are negligible. A 
correction to the upstream conditions due to vorticity effects down- 
stream could be made using the approximate theory suggested by 
Marble (Reference 15) for the three -dimensional case. If the hub 


ratio is large, the simpler two-dirnensional approximation given by 


Rannie (Reference 1) can be used. Once the upstream conditions have 
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been changed using this higher order correction, the corresponding 
new solution can then be found by the same procedure as previously 


outlined in Part III. 


B. Simplification of Flow Patterns 


An investigation has been made to determine whether or not 
there are any special cases in which the equations, used in this analysis 
of three-dimensional, compressibie flow, reduce to a simple expression 
for the axial velocity. Consider the integrated radial momentum 


equation: 


2 A 8 
E - 12 rA 2-4 22-8 20)42 «G (22), 


Since the value of P3 appears in tke integrated mass flow 
relation, Equation (21), to the power == » itis clear that, generally, 
no simplification can be made for a channel of arbitrary shape where 


the values of 


do 
O, 3 $, d 03 


d£, and dE cannot be expressed in 
terms of constants and the radius = : 

However, one simplification does occur in the example of a 
rotor in a cylindrical channel. In this case, if the velocity after the 
rotor is purely axial, Equation (22) reduces to P3 = CAT . NEP in 
addition, the total temperature and pressure are constant across the 
channel at station 2, and if Ne F A (vortex flow), Equation 
(16) for a straight channel reduces to $4 = Ci. ihe constant C, is 


evåluated by applying the energy equation at a boundary, and the mass 


flow integral becomes: 
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been changed using this higher order correction, the corresponding 
new solution can then be found by the same procedure as previously 


outlined in Part Ill. 


B. Simplification of Flow Patterns 

An investigation has been made to determine whether or not 
there are any special cases in which the equations, used in this analysis 
of three-dimensional, compressible flow, reduce to a simple expression 
for the axial velocity. Consider the integrated radial momentum 


equation: 
z v [qe 
- [Fra 
2 


Since the vaiue of P4 appears in the integrated mass flow 


2-3 de. y Hd e (22), 


relation, Equation (21), to the power P ; itis clear that, generally, 


no simplification can be made for a channel of arbitrary shape where 


ds d 93 
the values of a Z, and P; LE 


terms of constants and the radius = ; 








cannot be expressed in 


However, one simplification does occur in the example of a 
rotor in a cylindrical channel. In this case, if the velocity after the 
rotor is purely axial, Equation (22) reduces to Pr = CAT en in 
addition, the total temperature and pressure are constant across the 
channel at station 2, and if A,= A (vortex flow), Equation 
(16) £or a straight channel reduces to $4 * C. The constant C, is 


evaluated by applying the energy equation at a boundary, and the mass 


flow integral becomes: 
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or, since doe is constani, 
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Even in this simplified case, Equation (39) must be solved by 
trial and error for the constant C,. This particular example gives a 
constant mass flow per unit area and the work is taken out uniformly 
across the channel. 

With the exception of the above case, it can be generally stated 
that there are no special conditions or flow patterns in three-dimen- 
sional, compressible flow through turbomachinery using this method 


which result in a simple closed-form solution for the downstream axial 


velocity. 


C. Proposed Optimum Turbine Design Condition 

In the design of a turbine, the condition which normally controls 
tne amount of heat that can de released in the combustion chamber is the 
limiting temperature that the first stage rotor blade root can safely 
withstand in continuous operation. Since the blade stresses decrease 
from root to tip, the allowable temperature is higher at the tip than at 
the root. This fact allows the possibility of increasing the arnount of 
heat released in the combustion chamber, and hence the energy in the 
jet exhaust, for a fixed amount of work taken out of š turbine, 


In formulating the optimum turbine design condition, reference 
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will be made to a one-stage turbine with an arbitrary, but fixed, 
boundary configuration, such as is shown in Figure 5. Itis assumed 
that the total pressure in the combustion chamber, the blade speed, 
mass flow and total amount of work to be extracted by the turbine are 
fixed by conditions in the compressor. It is further assumed that the 
distribution of the maximum allowable total temperature across the 
channel at station Z, based on limiting blade temperatures, is known. 
Since the total enthalpy is conserved through the stator row, the limiting 
total temperature profile at station lis therefore known. 

in order that the maximum energy be left in the exhaust gases 
at station 3, the condition is imposed that the velocity shall be axial 


after the turbine ( A,= 0) and that the total enthalpy shall be constant 
dize 


3 
To simplify the formulation of the proposed optimum turbine 


across the channel at this station ( = 0). 

design condition, it is assumed that the known limiting total temperature 
profile at station 1 is linear, and that the combustion chamber can be 
constructed so as to give this desired total temperature variation across 
the channel. 

The design problem can then be formulated as follows: Givena 
one-stage turbine with fixed boundary configuration, combustion chamber 
total pressure, blade speed, mass flow, and linear total temperature 
profile at station 1, what are the blade spouting angles of stator and 
rotor so that a given amount of work can be extracted by the rotor and 
yet leave the total temperature constant and the velocity axial at station 
32 

With the above assumptions and conditions, Equation (27) can be 


written: 





sr, 


HG) = Sa àz (40), 


and Equation (10) in terms of total temperatures becomes: 


7 2 
2, 62r 322712 SE (41). 


If 7 2 is to be constant at station 3, then the work must be taken 
out across the channel in such a way that: 
A flee (42). 

Since the linear change in total temperature at station 2 is 
assumed known, the distribution of work across the channel is deter- 
mined by Equation (42) and the total amount of work to be extracted. 
The solution to this simplified inverse problem is found in a manner 
similar to that illustrated by Example X of Appendix F, and the blade 
trailing edge angles of stator and rotor are completely determined. 

In view of the discussion given in Part UI, Section B, on 
limiting blade angles, it will obviously not always be possible to obtain 
the above proposed optimum condition, particularly if the upstream 
total temperature gradient is large. Furtherrnore, the above consi- 
derations could not reascnably apply to a multi-stage turbine since 
boundary layer effects and mixing would secon destroy the original 
total temperature profile. However, the designer shouid not overlook 
the possibility of increasing the gross thrust of a turbine by the use of 
a total enthalpy gradient based on blade limiting temperatures. (See 


Example XI, Appendix F.) 
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V. RESULTS AND CONCLUSIONS 


An approximate method has been developed for the solution of 
three-dimensional, compressible flow in axial turbornachines. The 
methods and equations used are extremely flexible and can be used to 
solve the direct and inverse problems, subsonic or supersonic, isen- 
tropic or polytropic for either a turbine or a compressor. The upstream 
conditions, blade shape (as defined by the blade trailing edge angle), 
blade speed, and axial boundary configuration can be prescribed in the 
direct problem, while the upstream conditions, blade speed, blade 
loading, and axial boundary configuration can be prescribed in the 
inverse problem. In each case, six unknowns are solved for at the 
downstream station. These are the three cornponents of velocity, the 
pressure, the density, and the temperature. In the analysis, maximum 
utilization has been made of the given channel configuration. 

The equations are derived for a blade row with infinitely many 
blades, on the basis of isentropic (or polytropic), non-viscous, com- 
pressible, axially-symmetric, steady flow. Radial or nearly radial 
blades are assumed so that the radial force can be neglected. Separa- 
tion of flow and boundary layer effects are not considered in the analysis. 
For an arbitrary boundary shape, a special type of blading is assumed 
whose only effect on the flow is to turn it through a specified blade 
angle distribution across the channel. Using the assumption that the 
blade shape is specified by the blade trailing edge angle only, and 
baffling that this spouting angle is very nearly independent of condi- 
tions ahead of the blade row in the operating range, the analysis of the 


flow has been rade at the station immediately downstream of the blade 





-3- 


row in terms of specified conditions immediately upstream of the blade 
row. 

This method greatly simplifies the analysis of flow in turbc- 
machines, since, with the above assumptions, itis not necessary to 


solve for the complete three-dimensional flow field through the blade 





row, but only to consider stations immediately upstream and down- 
stream of the blading. 

The iteration process develcped is simple and straightforward, 
and the same general procedure is used in the solution of all suggested 
tyves of problezns. The calculations can be done with a slide rule, and 
the accuracy of the solution can be improved by increasing the number 
of racial steps across the channel. Corrections can be made for the 
displacement of the assumed streamlines and for the effect of down- 
stream — on upstream conditions to further improve the solution. 
It is believed that the results thus obtained by this approximate method 
are sufficiently accurate for design purposes. 

Once the boundary configuration and other necessary conditions 
have been assumed for any particular problem, experience in working 
various examples by the proposed iteration process has shown that a 
basic solution for the direct problem accurate to four significant 
figures can be obtained for a stator in approximately four hours and 
obtained for a rotor in approximately six hours, The calculations for 

—the inverse problem are somewhat simpler and less time is normally 
required for this solution. 

In order to evaluate the term = which appears in ooth 


the integral for the axial velocity and the integral for the pressure ia 
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the case of an arbitrary channel shape, it is necessary to estimate the 
slope and curvature of the streamlines at the station in question. It is 
to be noted that, for a reasonable channel shape, a relatively large error 
in this estimation will have little effect on the solution obtained. The 
term oe appears in the general axial velocity integral for the 
rotor (Equation (18) ) and for the stator (Equation (26) ) as a small 
correction term, and a large error in the assumed values of slope and 
curvature will have a small effect on the resulting axial velocity profile. 
This is especially true for the rotor where the term 2G is normally 
dominant. 

Considering the pressure integral (Equation (22) ), it is to be 


noted that here, again, the contribution of the term # is small 
2 





for a reasonable boundary shape. For the stator, the term is 
the predominant one, and the terms O ae and P ^ are 


usually small corrections. Because of the small magnitude of the factor 
multiplying the integrand of Equation (22) in the iteration, a relatively 
large error in E will have but a slight effect on the resulting 


pressure gradient. After the rotor, where the velocity is predominantly 


2 
axial, the term SAM might be of the same order of magnitude as 
do do s é 
= and = . As a result, it is possible to have a 
e az Pp ar p 


negative pressure gradient radially after a rotor if the slope and curva- 
ture of the channel is large. 

It is impossible to make any general statement as to the relative 
contribution of the integral (Equation (18) ) and the constant of integration 
Cy to the resulting axial velocity profile. The relative size of the 
integral and the constant may be of any magnitude depending on the 


channel shape, total temperature gradient, blade spouting angles, and 
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upstream ccnditions cf any particular example. 

Several interesting observations are possible from the results 
of the examples given in Appendix F. The channel shapes assun.ed for 
these iterated problems are given as Figures 3, 4, and 5. in each 
example, the assumptions, channel configuration, upstream conditions 
radial coordinates, and other data are specified, and the resulting 
solution at the downstream station is tabulated. Cne complete iteration 
is given in Exarnples Vb, Vliib, and IX b to show the method used in 
each case. 

Examples i and ll use the channel configuration shown in 
Figure 3. In both of these examples of the direct problem, the mass 
flow, blade speed, blade spouting angles, and total heat content in the 
flow are the same. Example I, however, has no radial total tempera- 
ture gradient at station 1, while Example Il has a large radial enthalpy 
gradient. The effect of a total temperature gradient across the channel 
in changing the axial velocity profile is clearly evident from the results 
obtained. Normally for a stator, it is to be expected that the positive 
pressure gradient would make the axial velocity higher at the root than 
at the tip. This is the case in Example la for no temperature gradient. 
However, Example IIa shows that the effect of the positive total 
gradient imposed is to make the axial velocity higher at the tip than at 
the root, and reduce the Mach number slightly at all stations across the 
channel. Also, the resulting total pressure at station 3 has been in- 
creased by the use of the radial enthalpy gradient. 

Examples IV and V use the channel configuration of Figure 5, 


and have the same mass flow, blade speed, blade spouting angles, and 
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conditions at station 1l, Example IVa is the subsonic stator solution 
while Example Va is the supersonic solution for the same stator. 
Examples IVb and Vb are the subsonic and supersonic rotor solutions, 
respectively, for conditions after the subsonic stator of Example IVa. 
The mass flow functions for the subsonic and supersonic rotor examples 
are plotted in Figure 6, From this graph, the position of the displaced 
streamlines has been estimated and plotted on Figure 7, This plot of 
the streamline shift in each case, clearly indicates that the shift is in 
one direction when the "governing" Mach number is subsonic, and in the 
opposite direction when the "governing" Mach number is supersonic. 
The "governing" Mach number has been defined by Monroe (Reference 7, 
page 51) as the Mach number based on the meridional velocity instead 
of the total velocity. The results shown in Figure 7 clearly agree with 
Monroe's observation that the deflection of the “meridional Mach 
surfaces" is in different directions for subsonic and supersonic 
"governing" Mach numbers, and that the deflection of the streamsurfaces 
will be less if the governing velocity is transonic than if it is entirely 
subsonic or entirely supersoníc. 

Examples VI, Vil, and VIII use the channel configuration of 
Figure 5, the same conditions at station 1, and the same blade speed, 
blade spouting angles, and mass flow. Example VI is the isentropic 
case, Example VII is a polytropic solution with n= 1.37 (7 = 94,5 ings - 
cent), and Example VIII is a polytropic solution with n - 1.33 
(7 = 86.8 percent), From these examples, the effect of a decrease 
in "polytronic efficiency" on the three-dimensional velocity field and 


distribution of energy in the flow can clearly be seen. When the 
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"efficiency" has dropped to 94, 5 percent, the average axial velocity at 
station 3 has increased 2.7 percent, the average total pressure has 
decreased 2.5 percent, and the work output has increased 4.7 percent 
when compared with the isentropic solution. As the "polytropic 
efficiency" drops to 86.8 percent, the results of Example VIII show that 
the average axial velocity at station 3 has increased &.2 percent, the 
average total pressure has decreased 8.0 percent, and the work output 
has increased 17.0 percent when compared with the isentropic solution. 

The results of these examples clearly show the reason for the 
"loading up" of the first stages cf a multi-stage turbine as the efficiency 
drops. As the efficiency decreases, more work is taken out in the 
initial stage, and less total pressure is left for the remaining turbine 
stages. Thus, the last stage of a multi-stage turbine may well act as 
a compressor if the efficiency of the system has dropped sufficiently. 

In view of this phenomenon, it is obvious that the design engineer should 
design his machine for the expected efficiency in the operating range, 
and not on the basis of isentropic flow. 

Example XI has been constructed to illustrate the increase in 
gross thrust obtainable by the use c£ a radial total temperature gradient. 
Assuming a one-stage turbine configuration similar to that shown in 
Figure 3, and keeping the mass flow, combustion chamber total pressure, 
blade speed, and work extracted constant in each case, the gross thrust 
of the jet exhaust has been calculated with and without a total tempera- 
ture gradient. In order to give the example more physical meaning, 
blades were used which were tapered in cross-sectional area from hub 


to tip. Having calculated the centrifugal stress at various radial points 
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along the blades, and using the high temperature data for ''Vitallium" 
(1.5. 21) (NR-10) given in References 16 and 17, it was estimated 
that the rotor blade root could safely stand a temperature of 1450°F, 
while the blade tip temperature could be 1850%r, 

Using this data, the gross thrust oí the turbine was calculated 
for: (a) A constant radial total temperature of 1450%F, and (b) a linear 
total temperature gradient of 400%F with the same root total tempera- 
ture of 1450°F. The results of these examples show that an increase 
in gross thrust of 917.8 pounds, or 7.7 percent, was obtained by the 


use of the radial total temperature gradient indicated. 
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APPENDIX A 


SYMBOLS AND DEFINITIONS 





Be - Stator spouting angle measured from the plane perpendicular 
to the axis. 
Ba - Rotor spouting angle measured relative to the rotating whecl. 
€ - Axial distance between stations in feet. 
Ci - Constant obtained from axial veiocity integral. 
C, - Constant obtained from radial momentum integral. 
Cp - Specific heat at constant pressure. 
A - Used to denote incremental changes of a quantity in the radial 


direction. 


F( Es) - Cotangent of stator spouting angle Pe as a function of radius. 


ES ame 


G (53) - Cotangent of rotor dier n Spouting angle Ps as a function 
of radius.  G/§,) = a Matiz 


Y - Ratio of specific heats. 
g - Acceleration of gravity (32.17 ft/sec”). 


H(5, - Dimensionless work output of rotor as a function of radius. 


| H(£,) * $202 * $325 











i Wu 
m - Dimensionless mass flow. m= 2179 Py Re 
M - Mach ber based on total velocity. = 
ach number based on total velo me! = ee 
Mg - "Governing" Mach number as defined on p. 87 of Reference 7. 
M _ Dc ee 
» rT 
n - Polytropic exponent. 
ee ca 24 _ ZN 
7) - Polytropic efficiency. h e 
U) - Angular velocity of rotor in radians per second. 
p - Pressure in pounds per square foot. 
22 


P - Dimensionless pressure. P = -z 








J un N 


«X 


W 


3% 





- Density (lb. sec. 2 lft. 3. 





. t | «?ASe 
- Dimensionless density. = EN 
2 
- Gas constant (1715 af - 
sec^?r 


- Reference radius in feet. 


- Radius from the axis in feet. 


è A z Za 
- Dimensionless radius. 5 EE 





- Absolute temperature in degrees Rankine. 
AR 
- Dimensionless temperature. C = GERZ 7 


- Radial component of absolute velocity in feet per second. 


U 


- Dimensionless radial component of absolute velocity. 0= AS 


- Tangential component oí absolute velocity in feet per second. 
- Dimensionless tangential component of absolute velocity. A = Ser 

o 
- Axial component of absolute velocity in feet per second. 


- Dimensionless axial component of absolute velocity. ¢= a 





- Dimensionless total absolute velocity. 8 = \|¢*+ A°+ 0° 

- Mass flow in pounds per second. 

- Axial distance along axis of rotor from reference point in feet. 
- Dimensionless axial distance. ¢ = 5 

- Subscript to denote station immediately - of stator. 


- Subscript to denote station immediately downstream of stator 
and immediately upstream oí rotor. 


- Suoscript to denote station immediately after rotor. 
- Subscript to denote inner boundary of channel. 

- Subscript to denote outer boundary of channel. 

- Slope of streamline. 


- Rate of change of slope of streamline. 
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APPENDIX B 
Note: In Appendices B, C, D, E and F, Equations of lower number 

than 43 refer to relations previously given in the text. The 

numbering of equations is consistent throughout the thesis. 
Development of Zquations for Direct Problem (Isentropic). 

The equations are derived for a rotor with infinitely many blades 
on the basis of isentropic, non-viscous, compressible, axially-symmet- 
ric, steady fiow. Radial or nearly radial blades are assumed so that 
the radial force can be neglected. The channel shape, blade shape, 
mass flow, and upstream conditions are assumed specified in any par- 
ticular example. 

Tine derivation follows the general procedure given in Refer- 
ence 1, Section 2:2. 2 ior the incompressible case. The energy equa- 
tion is written for any two points (for example, b, and bz» Fig. 5) on 
a streamline suriace between stations 2 and 3. Using absolute veicc- 
ities, the equation in dimensionless form is: 

2 
= $2 A x E E "(at Ef ag |- z $3 (43) 


where 


2 
yg Ter = 77 La (0 t Aet $; ), 


and 


/ E 2 2 
7 = AT + —/G, + + 
7-, 6st Fp ts + BO Ay 9; ) 


The energy equation is now differentiated in respect to the rad- 


ius at stations 2 and 3: 





pa 
Ct 





* dist d 
7- d£, df 


£z (f Jie = 1255 lar) zz EA), ag, 9 


Using the gas equation, 





A 

s ES , the first term of the right 

side of Equation (44) becomes, with the isentropic relation 
-L M 

a , 

Ya Pal \% -— 
zdl _ zd 
rı di — 


a [Al _ 7 q Pa? 7) 
Y=/ (2) - 7-1 2 (2 ^ 


Using the differential form of the continuity equation for a 
Stream cylinder oí small radial extent 


Po Qa p d 3 


(45). 


$5 4, 5, a £, 


(8) 
and difííerentiating Equation (45) as a product, the following is obtained 


Ir 
S 
— 

D Io 
Sa 
! 
x 
Nik 
a 
qo 








2 d hj 
"A 
3 77% dé, 
4 


asco bs ds §3 (15) 
dé, Pe És 
Using the radial momentum equation at station 3, 
AA „20.2022 
45 26; 53 





E g 803 
2 E zuo $3 
and combining 


(11) 
quations (8), (11), (15) and (44), the following differ- 
ential equation results: 





46 
LE re a v) ee, ET 2, gn). ge 
7- Br de, Paño 2 df E 7) d é; 3 


= 7. cte _ d($2A2) 
i / le, d É las, 





P AGE? EA 32" EICHE, i 


AY ded or 3e 2X) v Qe A 2 (2) (46). 
Pa! Pea [71 DE 98, FI Q5 Sm Qe 


For the direct problem it is assumed that the rotor spouting 


is known as a function of radius at station 3. It is there- 


angle P3 
fore possible to describe the relation between A; and $3 along 


the radius with a function G( $3) such that: 


Az = 6 ($,) $,- É, 


(17). 


Eliminating Az from Equation (46) with this relation the fol- 


lowing differential equation is obtained: 


Gt- (64-6), E (ee, £) ej ag, BGs £) 


"% e G 
IE ta e E) nre 





G' os *6 3 (e 9) * P 


z 
(es 2 E Z oe 9) o as. (. Jz 28, 26 


Pesa 7219 fa ? fz d 











OP; tg ee = a(1rG2) 
Y £e) 2(6%) 2f, 


93 
le 
“(A ss Er, N, sE t reg 647). 
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803 _ à 
EP, 3 > SE and the known quantities at station 2 are 
considered to be functions of E 3 (to be determined in the iteration 


process) the integral for Pz can be immediately written down as: 





5 d $ £ 

+f GF 3 en. = 

bor £^ G^) L Bie y (2) É 2005. | d ($e 2) 
p= $, 733 S AG? fe] Pee |7 d f> ET 





ein 


EP ff 
- Zu) 264) + E [26+ £2 26) e d£ 4C 


(18). 


An adternate form for the square bracket in the above inte- 
grand can be obtained by using the energy and radial momentum equa- 


tions at station 2: 


in the case that Pt = 1 (isentropic expansion to station 2 from 
stagnation conditions), the bracketed term in the integrand may be 


further simplified to: 





7 die zZ) -^ 
z a, (17% 7) - $e (fade) (20), 
The integral form of the mass flow equation is: 
^ Ps 7 
n -S eR] ahas, (21), 
¿ 


and the integral form oí the radial momentum equation is: 


; 5 
= TAE de Az _ des. 7 
— "F CA ra Edere (22), 
5 


Equations (10), (13), (17), (18), (21), (22), and (24) are used 
in the solution of the direct problem for the rotor (isentropic). 


In the case of a stator, Equation (17) becomes: 


F(&,) bz (25), 


G is replaced by F in Equation (18) and the terms zh Az) and 

2 
2G are omitted. Equation (18) becomes Equation (26) for the isentropic 
stator. Equations (13), (21), (22), and (24) (with appropriate sub- 


scripts) remain unchanged, and the energy relation becomes Equation 


(9). 
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APPENDIX C 
Development of Equations for Direct Problem (Polytropic). 
The equations for a polytropic process through a rotor are 
derived with the same assumptions as used in Appendix B with the 


exception that the isentropic relation (13) is now replaced by: 


&- (afte à 


where n is related to the so-called polytropic efficiency 7 íor an 


expansion (n < y) by: (see Reference 12, p. 449) 


n A 
an S (49). 


The energy equation (10), the continuity equation (8), the rad- 
ialmomentum equation (11), and the blade relation (17) remain the 
Same in differential form. 


The first term of the right side of Equation (44) now becomes: 





But, 
a 445 zr A dQ@s _ Lord Es 
094 Qadi | 74 Q^ a£, a OS 
So, as before, 
? d a) _ L 4l, pd (E (50) 
NES O 713 ALE \ @ l 
-l a£ Q; q, d 3 3 I 
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The right side of the differential equation is obtained in the 


same form as Equation (46) with a correctioa term as follows 








AF _ (er) 22 
= 77 
Q3 Q3 ) 2 
4 
x pd [R L pr pad 
x 7I (oe = 774 É á 
774 7 ae \Q,/ = 7-79 E dé, 4). 4 % igh (51) 


2 


Using the differential form of the continuity equation (8), Equation 
(51) finally becomes: 


= 
7 





= L a z 22-4 
PAR ]- ZA (BVA poh em (BP) a 
dí, 5 7/3 LEN Ge Pal Pa fo (M) 


The differential equation corresponding to Equation (46) for 
the x ee case now is: 


ag 9S, 


e 
- 3 Do E Y ler 2($5A 2) 
B. * $3: * )* 3; Tm Tee 


7 


2-1 
7 pi O 96 (9-2) zu 
7; ^ 2£, 4) a er ? ($25 )91-.) 2027) Ge ) 2A 75 (53). 


The relation corresponding to Equation (47) becomes 


295 +g on ll alt?) _ af 9% J_r 3%: _ 3629.) 
2£, 3 lEll6)  2(1+69) 95, - 


Pa fa | 86. 96, 
-/ 25 22-/ 
nod x s = 265 7-2 22) 28 ET 
ae 38, )- 26, RN SS TH (54), 


The integral equation corresponding to Equation (18) is 





5] 











fs é 
d£ 3 pfe 
mi 3 f 50 i 
i z 5 8)| | dicen ; (ay £ ra dtere A(Ge Az) 
£, N/^6? (^6? f | pafe | o! Ll $e 


2-1 


Y pra par / do; 1 en) (02)\dB dé 
EET: E) leo de * $i) ali d y A 





(55). 


If Pre = l at station 2, the square bracket under the integral 


will be: 


7 Gee L m) 
E dé, (^-^ / d£ | (56), 


The integrated form of the mass flow equation is 


fo 4 
-S a(R] as (57) 
a 


and the integral form of the radial momentum equation is: 


72 


5 72-1 
- e 3 ga 2% a 685 
get = (e E “2, "ag 45 C. -— 
£. 


Equations (10), (17), (24), (48), (55), (57) and (58) are used 





in the sclution of the polytropic direct problem for the rotor. 

In the case of the stator, Equation (25) is gen in place of 
Equation (17), and the terms iz. Az) and 2G are omitted in 
Equation (55). Equations (10), (24), (48), (57) and (58) remain un- 


changed (with appropriate subscripts). 
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APPENDIX D 
Development of Equations for the Inverse Problem (Isentropic). 
The following equations are derived with the same assump- 
tions as given in Appendix B with the exception that now the blade 
relation, Equation (17), is replaced by the following given function 
which describes the work output of the rotor as a function of the 


radius £; : 
E SDR (27). 


For the isentropic case, the derivation remains unchanged 


up to an including Equation (46), which may be written as: 


w fv 


À LU = re AMT 7 el 
1 ls) AGA ZE 


3 
MA E 
7,3 gg iq 39 €. 
fa 7 (28). 


Using Equation (27) to eliminate A 3, the differential equa- 


tion corresponding to (47) is developed as follows: 


vn 


a, 1dAs_ Led dlEar,) 


/ 
3 eds, F £ zi. 


a Ps JAs d es P. 3753 d 72 a EA. 
az (+40) Latin) - (af tos Sp 2o 


e 











d _ dH_ (BY bh d 
2. (H-£ Aa) ^. dé, E be £, 149292) (59) 





A 








123 
+(e EEA iz E ut AZ) 


u \ 
ds _ (E Aus Eu C NC d Bre 
d£, E 


N 
A 
S 
N 
\ 
A 


Assuming that the right side of the above Equation (60) is 
a known function of 5 3 (to be determined in the iteration process), 


the integra! equation for the axial velocity can be written as: 
fs 


PNTE d. a 222) dl 
$7 [255 Fi af we Z^ Ce 4o) tE 


z y Ad 27 af, 








4 


(60). 


C3 


(29). 


Since the mass flow integral and the radial momentum inte- 
gral remain the same, Equations (10), (13), (21), (22), (24), (27) and 
(29) are used in the solution of the isentropic inverse problem for 


the rotor. 
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APPENDIX E 


Development of Equations for the Inverse Problem (Polytropic). 


The following derivation is made assuming the same conditions 
tor a rotor as given in Appendix C. The work output of the rotor as a 


function of radius is again piven by Equation (27). 


The differential equation (53) still holds, which is: 











2 
As II A 2 E z 9720 — 2$222) 
"MEE TUR ME DRE VES: P S 
-— ^o 2980, 9-22 2= 
7-1 4 IE, Ass A TT Qə le ) (53). 


The resulting integral equation is similar to Equation (29) with 


the addition f a correction term: 





2) xdle + paid Wied ge Àz LO; 
E E, 21 de 2» 2 (Aj ) rc aay “dE 
i 2 
(7-72) d/27 + HE are ie 
P; a (2 dm [1+ zn z pdg C (61), 


The mass flow and radial momentum equations remain the 
same a8 in Appendix C, and Equations (10), (24), (27), (48), (57), 


(58) and (61) are used to solve the polytropic inverse problem for 


the rotor. 





AF PENDIA F 


Iterated Examples 
The following data are used for Examples I-X inclusive: 


lb 


W = 734 rad/sec (7006 RPM) p,, = 16,300 — Rz 1.29009 £t. 
Lt (t O 
vo = | 240 wR_ = 917.5 ft/sec C = .2396 — nU 
3 S p bor 
R =1715 E - ia 
sec” F en? = ‚8416 x 10 ——. c= .250 ft 
sec 


Example la. 


Solution for the stator of Fig. 3 (cylindrical channel), with no 
radial total temperature distribution and coordinates for equal areas 
in the plane perpendicular to the axis. This is a solution for the direct 
problem (F (15 ) given) for the isentropic, subsonic case. 

For the stator with F constant, P,, - l, and 0 - 0, Equations (18), 


(21), (22), (13), (25) and (43) of Appendix B become: 


R 











fe E 
7 Ll 
e- o. - 16) + 2 ag (I-A P) £ real y se (62) 


/ E J+F? 


3 
N 


^s 4 
[aly esas, (63) 





(64) 


NO 
i 
FF el 
N | T 
N 
Ir 
SN 
eI. 
R 
Ro 
+ 
C 
zl 


$ - (Af'- (&* d 


À = F $ (F = const.) (66) 
yy = ED «i (Xe) (67) 
y Ma ^  4*e 2 2 2 
M = $ (Mach number) (68). 
Yzz 
Conditions at station !: (given) 
yt. 00 ©), = 23028 l,,= 4.1968 $, = .4857 
F, = .9723 C, = .23354 l, = 4.1631 A,= 0 
E 0 F=2.75 m = „020417 we er, - 2er 
£ sec it 
Solution at station 2: 
a, b, C5 d, e, 
a 8000 ‚8544 .9055 .9539 1.0000 
$5 .7926 7479 7105 6785 .6508 
Aa 2.1796 2.0567 1.9539 1.8659 1.7897 
P 4928 .5365 ‚5730 ‚6041 .6307 
Tə 3.4283 3.5128 3.5794 3.6340 3.6789 
P 14372 ‚15274 ‚16007 ‚16623 „17144 
M 1.059 .986 .929 .880 .841 
AC,  *.0001 - .000Z -.0001 -.0004 - .9002 
ae RES Poe Fag Gil Cg Og 02 20 
C, = .8169 


57 
Note: The column headed A Z 2 above represents the difference 
between T > as calculated by the iteration process and lo 
obtained by applying the energy equation along the stream- 
lines indicated after the solution has been reached. The same 
applies in the solution of the following examples. This is a 
rough check on the accuracy of the solution. Obviously, if an 
exact solution had been found, AZ would be zero. 
Example Ib. 
Solution for the rotor of Fig. 3 (cylindrical channel) with no 
radial total temperature gradient and radial coordinates for equal areas. 
This is an isentropic, subsonic solution for the direct problem (G { LE ) 
given) for a rotor. [ 
For the rotor with G constant and Pot = l, the applicable 
relations are Equations (63), (64), (65) (with appropriate subscripts) 


and the following three equations: 








£ PG? 
3 


s í 
/ P 7 = Em C, 
a AA Pte] seo gag |o s (o 
/1 G*) 
$ 


^; a C%- $, (G = const.) (70) 


HTM AE (AA) A, (71). 


Conditions at station 2 are the same as after the stator, Example la. 


G= 1.411 m = .020417 


Solution at station 3: 


“3 
b. . 8000 
$5 . 6624 
Az . 1346 
P. . 5862 
7, 3.6025 
Q3 . 16269 
Pa «6244 
Tj, 3.6679 
Qa . 17019 
AL, 0 
Cira ¿5710 


1 
Example Ila. 


Ba 
8544 
‚6779 
‚1021 
. 5865 

3. 6034 
‚16278 
6254 

3. 6699 
‚17041 


-, 0007 
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= 

. 9055 
. 6955 
. 0759 
. 5865 
3.6034 
. 162176 
. 6264 
3.6717 
. 17061 
-.0017 


C, = .8585 


2 


ds 
‚9539 
‚2146 
‚0544 
, 5865 
3.6035 
16276 
6273 
3.6734 
‚17077 
-, 0035 


x 
1. 0000 
. 7346 
. 0365 
. 5865 
3.6035 
. 16277 
. 6284 
3.6751 
, 17099 
-.0058 


Solution for the stator of Fig. 3 (cylindrical channel) with the 


same conditions and assurnptions as given in Problem la except with a 


radial total temperature gradient. 


Equations (62) - (68) of Example la are applicable. 


Conditions at station 1: (given) 


al 
f, . 8000 
Te — 3.3819 
E. 1000 
Qi, . 29569 
P, - 9723 
T, 3.3548 
Qi . 28981 
$, 4355 
A, 0 

F 38 2.75 


bi 


8544 


3.8252 
1. 00 

. 26142 
. 9723 


3.7946 
. 25622 


. 4628 
0 


» 
. 9055 


4. 2415 
1.00 
A317 
.9123 


4.2076 
. 423108 


. 4871 
0 


= „020417 


di 


. 9989 


4.6361 
1.00 
.21570 
9123 

4. 5990 
. 21141 
. 5096 
0 


a Tip 
a 


= +8. 1495 


"M 
1. 000 


5. 0118 
1.00 
- 19953 
ES 
4.9717 
. 19556 
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Solution at station 2: 


as Ba Co da ez 
5s . $000 . 8544 .9055 9539 1.0000 
Pa .6778 . 6805 . 6829 ‚6848 . 6864 
Aa 1.8639 1.8714 1. 8780 1. 3832 1.8876 
P, . 5295 . 5713 .6061 . 6355 . 6611 
Z > 2.8200 3.2599 3.6762 4.0729 4.4527 
Q, . 18777 . 17525 . 16488 ` . 15602 . 14845 
M . 999 . 930 . 883 . 842 .805 
At, 0 -.0012 -.0052 -.0104 -.0172 

C, = .5565 C, = . 8339 
E? TOL, Q25 Aje. 9279 


Exarnple lib 
Solution for the rotor of Fig. 3 (cylindrical channel) with the 


—same conditions and assumptions given in Problem Ib except with a 
radial total temperature gradient. Equations (63), (64), (65), (69), (70) 
and (71) of Example I (with appropriate subscripts) are applicable. Condi- 


tions at station 2 are the same as after the stator, Example Ila, above. 


G s 1.411, dlet 2 8. 1495, m s ,020417 
Solution at station 3: s 
2i b3 dE ds e3 
E . 8000 . 8544 . 9055 . 9539 1.0000 
$> . 5377 .6174 . 6920 . 7628 . 8303 
Aa 7.0413 . 0168 . 0709 . 1224 . 1716 
P, . 6009 . 6009 . 6009 . 6009 .6012 
T. 2.9397 3.3072 3.6670 4,0081 4.3334 
Q, . 20551 . 18168 . 16386 . 14990 . 13871 
Au 0313 , . 6381 . 6422 . 6447 . 6461 
Q4,  -21289 . 18964 . 17184 . 15759 . 14603 
tz, 4.9655 3. 3643 3. 7373 4. 0895 4. 4235 
aT, 0 +.0035 +.0012 -.0039 -.0126 
C, 2.4635 C, = . 8646 


1% 2 
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. 
Aíter correcting the upstream conditions at station 2 for the 
displacement of the streamlines from the position originally assumed, 
the following corrected solution at station 3 is obtained: 


(See Problem 


Vb for method. ) 


az b3 E d, 83 
M 8000 . 8544 .9055 .9539 . 0000 
$3 .5311 . 6106 . 6853 . 7560 . 8234 
Az . 0506 0072 0615 1128 . 1618 
P, . 6029 . 6029 . 6029 . 6029 . 6032 
T; .9268 3.2771 . 6296 3.9868 i. 3374 
Q3 . 20602 . 18411 . 16623 .15127 . 13917 
E, -6827 . 6373 .6417 . 6453 . 6476 
EL 2.9675 3. 3303 . 6951 4. 0648 . 4263 
Oy, «21325 . 19166 . 17381 15879 . 14641 
AT; 0 -. 0001 -. 0022 - 0055 -. 0118 
C, 2.4578 C, = . 8654 


Example lila. 


Solution for the stator of Fig. 4 (divergent channel with straight 


walls} with a radial total temperature gradient. This is an isentropic, 
subsonic solution for the direct problem for a stator (F (5, ) given). 
For the stator with F constant and P,, = 1, Equations (63), (65) 


and (56) are applicable in addition to the following four relations: 





oy, 
+ 


GQ 
LEES (12) 


1+F 
2 


d = | EZ x E = 27)» a0 PS 
i ge tal li ) a, [7 








2 
% 7-1 
= 2-/ / Ale a e. $ a. 
Fra) lat tajo 
e 
D^ d / e 2 2 
ES le = T,“ + z (& + A: + P, 
202 . p eE 
24 IFNIEL “Nake 
Conditions at station 1: (given) 
a] m e dd 
5, .80 .85 . 90 .95 
T. 3.9931 4. 0950 4. 1968 4.2987 
Ay . 25043 , 24420 . 23828 . 23263 
Pie 1, 00 1.00 1. 00 1.00 
Py . 9696 . 9696 . 9696 . 9696 
T, 3.9581 4.0591 4. 1600 4.2610 
eu . 244971 . 23888 . 23309 . 221560 
6j , 4944 , 5014 ‚5069 ‚5112 
9, -, 0247 .0000 . 0253 .0511 
A, 0 0 0 0 
F = 2.75 e v^1.0378 5.028639 we 
/ 


i 
138. 4 


(73) 


(74) 


(75). 


el 
1.00 
4. 4006 
22724 
1.00 
. 9696 
4. 3620 
. 22229 
‚5142 
‚0771 
0 


lb 


li ls 


Scc 





Solution at station 2: 


a 


b 





2 2 3 en e 

E. . 792 . 850 .910 970 .30 
$5 . 1410 . 7046 . 6135 6469 . 6236 
Aa . 0377 9376 8521 1.7790 7149 
Oo . 0370 0 . 0337 0647 0935 
E. . 5247 . 5716 . 6122 6472 6711 
a . 3212 . 4904 . 6479 3. 7966 23376 
a, . 15798 . 16378 . 16782 . 17051 . 17210 
M 1.004 .930 . 876 .817 . 7179 
B . o -.0027 -.0061 -.0103 -. 0138 

C} =.6017 C,=.8317 Py. = Pi. O79, &=Z, 


Example Ib. 


Solution for the rotor of Fig. 4 (divergent channel with straight 
walls) with a radial total temperature gradient. This is an isentropic, 
subsonic solution for the direct problem for the rotor (G (£,) given). 

Equations (63), (65), (70), (73), (75) (with appropriate sub- 


scripts) and the following two relations are applicable: 





bs 
gb: 2 
B e £5 r alae 7), d ) We} one, URC 
= —p— | | 4/52) -5 LE ee zr d PU + “Ba 
+ g otyg?) 16 $. £7 d£ A 4g ev + dé 5 ad ¿ño (76) 
3 


xX 7 / (77). 
lache zart +2(8 Ae), 


Conditions at station 2 are the same as after the stator, Example 


Ula above. 
dTor 
Ls, 





G=1.41l m = . 024639 = + LEJT 8 Zee 
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Solution at station 3: 


a, 55 C3 d3 ©, 
6 . 780 . 850 .920 .990 1. 060 
P3 . 4567 . 4987 . 5423 . 5873 . 6347 
03 -. 0228 0 . 0271 . 0587 . 0952 
À2'  -.1356 -.1463 -.1548 -.1613 -. 1644 
P3 . 6500 . 6500 . 6505 . 6508 . 6508 
Q; . 18409 . 17954 17525 .17119 .16720 
T3 3.5308 3. 6209 3.2117 3.8027 3.8923 
Ps ‚6712 .6749 . 6780 . 6804 . 6828 
Q3 . 18836 . 18442 . 18051 . 17672 . 17304 
Tap 3.5634 3. 6600 3.7560 3.8514 3.9459 
A7,  *.0001 +. 0005 -. 0013 -. 0048 -. 0092 
C., 3.8871 C, = .8842 


Example IVa 


Solution for the stator of Fig. 5 (divergent curved channel walls) 
with a radial total temperature gradient. This is an isentropic, sub- 
sonic solution for the direct problem for a stator (F (§,) given). The 


same equations as used in Example llla are applicable. 


Conditions at station 1: (given) 


a) bi Ci di e, fi 

d . 80 .84 . 88 .92 . 96 1. 00 
7,, 83.7894 3.9524 *.1152 4.2783 4.4413 4. 6043 
E 1.00 1.00 1.00 1.00 1.00 1.00 
Q,, "26389 .25301 .24299 .23374 .22516 . 21719 
m .9400 .9400 .9400 .9400 . 9400 . 9400 
T, 3.7231 3.8832 4.0433 4.2034 4.3636 4. 5237 
€i .25248 .24207 .23248 .22363 .21542 . 20780 
$, .6804 .6960  .7090 . 7205 . 7293 . 7365 
Ə} -. 0340 0 . 0354 . 0720 . 1094 . 1473 
A 0 0 0 0 0 0 

F =2.75 atte 2.0786 m -.029242  W - 164.2 1P. 


ds sec 
i ¢ = . 260 ft. 





Solution at station 2: 


a, b, C5 d, e, t5 
do .785 ‚842 899 ‚956 1.013 1.070 
e, 47051 .6763 .6531 .6344 .6192 .6072 
A2 1.9390 1.8598 1.7960 1.7446 1.7028 1.6698 
9, -.0705 .0135 .0914 ‚1649 ‚2353 .3036 
P, .5417 .5860 6222 ‚6518 ‚6761 .6962 
T, 3.1805 3.3928 3.5937 3.7860 3.9713 4.1519 
G, .17022 Ames ‚17313 ‚17217 .17022 ‚16767 
M .977 .909 .854 ‚809 ‚776 ‚747 
Ale 0 0 -.0013 -.0039 -.0069 -.0118 
C, = .5694 C, = .8393 
T pt lp * lox z Lie i Car - S 


Example [Vb 


Solution for the rotor of Fig. 5 (divergent curved channel walls) 
with a radial total temperature gradient. This is an isentropic, 
subsonic solution of the direct problem for the rotor (G (5, ) given). 
The same equations as used in Example IIIb are applicable. 

Conditions at station 2 are the same as after the stator, " 


Zxample IVa. 


d lət 


G = 1.730 = 2.8593 m = .029242 
dg 7i 


Solution at station 3; 
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e 


23 ^ 
53 . 760 . 848 
bs . 1958 . 2641 
À3 ~. 4213 -. 3911 : 
05 -. 0294 . 0106 
P, . 6948 . 6981 
T, 3.4149 3. 5669 
Q4 , 20346 . 19573 
e 3. 4660 3.5998 
ES 21117 . 20027 
PES - 7220 . 7209 
M .212 .210 
al, 0 +. 0010 : 
C, =.1594 
Example Va. 


3 


. 024 

. JOSE 
. 3576 
. 1617 
. 4027 
. 8652 
. 18167 
2. 
. 18619 
. 7274 
. 236 

. 0049 


9064 


. 9012 


l 


4. 


$3 


112 1. 


. 4458 


. 2119 
. 7010 


. 17466 


. 17950 


. 7283 


- 266 


0568 4. 
. 17272 
< LaSi 
FOZ 

.0118 -. 


200 


. 29117 
. 3408 -. 


3148 


. 4024 
. 6927 
.0126 4. 


1461 
16707 
2018 


0196 





Supersonic isentropic solution for the stator of Example IVa 
(divergent curved channel walls and radial total temperature gradient). 
The conditions at station 1, the data, and the applicable equations are 


the same as for Example IVa. 


Solution at station 2: 


a, b, E co d, e, f, 
$5 . 9303 . 8906 . 8601 . 8371 . 8201 , 8087 
^a 2.5583 2.4491 2. 3653 2.3020 2.2553 2. 2239 
ĝ2  -.0930 . 0178 . 1204 . 2176 . 3116 . 4043 
P, 3171 3715 . 4169 . 4547 4856 .5102 
Ta 2.7294 2.9784 3. 2051 3.4157 3.6131 3.7990 
Gp . 11617 , 12471 . 13007 .13310 .13440 . 13430 
| M 1. 367 1.274 1.188 1.123 1.080 1. 038 
Als 0 +. 0039 +, 0034 -.0012  -.0082 -.0178 
C,-7.751à (6, 54909 ME HP 


2t = iy lo, z e, 


1 t 


2t it’ 
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=xarnple Vb. 


Supersonic isentropic solution for the rotor cf Example IVb 
(divergent curved channel walls and radial total temperature gradient). 
The conditions at station 2 are assumed the same as after the subsonic 
stator of Example IVa. The data and equations used are the same as that 
used for Example IVb. 

One complete approximation for this problem is given below to 
show the method of iteration. The procedure of correcting for the dis- 
placernent of the assumed streamlines is then given to show the method 
for this correction. 

Ihe iteration is started by assuming an axial velocity distribu- 
tion across the chanael. The tangential velocity is then calculated using 


Equation (70). The value of jv is found from Equation (75) and 
3 





the pressure distribution from Equation (73) and the application of Equa- 
tion (77) at the inner channel boundary. The axial velocity distribution 
is then determined using Equation (76) and the constant is determined by 
the mass flow Equation (63). The iteration is continued until the values 
of $4 assumed agree with those resulting from the iteration process. It 
is important to note that in the supersonic case, the iteration diverges 
rapidly away from the solution. Thus if an axial velocity distribution is 
assumed that is greater than the solution, Successive iterations will in- 
crease E until the pressure F3 has reached zero. If an axial velocity 
distribution is assumed which is less than the solution, successive 
iterations will result in decreasing values of $5 and eventually the subsonic 
solution will be reached. The supersonic solution is found by a simple 


inter polation once the solution has been "straddled". 
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The reverse is true for the subsonic solution for a given setup. 


The iteration process will converge toward the solution from any 


reasonable assumed distribution of $3: 


Data used in the solution: 








A4, = .088 Fr Ege = 10.0075 G=1.730 2G = 3.460 
2 
E 
=.7496 
"m 
a3 b3 C3 ds e3 f3 
.760 . 848 .936 1. 024 1.112 1.200 
2 -.150 .040  .230 ‚420 ‚610 — .800 
27/3 
Bal 
2E -.250 ‚100 „450 ‚800 1.150 1.500 
24 7/3 
Factor ‚02117 ‚02125 .02133  .02138 .02144 
2 ' 1.4310 1.5415 1.6411 1.7305 1.8098 
bE, 
/ 
Zaz 1. 2283 1.1316 1.0509 .9824 .9235  .8722 
3 
= (E$ 2 
Q, 6f) ‚01299 .01357 .01489 .01694 .01719 
d = ww = = EI 
- plade) . 7702 .8526 -.9333 -1.0018 -1.0842 
7- 1 u 
FE lb £,)  -006498 - 005993 .005768 .005560 


("Factor" above = 


gfe me (146°) 


Try C} = 1.3920 
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ag " ES 
Pz 1.7098 1.6727 1.6943 
Gh, 2.9580 2.8938 2.9311 
As . 2.1980 2.0458 1.9951 
O, -. 2565 . 0669 . 3897 
2 
As 4.8312 4.1853 3.9804 
5 
2 6.3568 4.9354 4.2525 
5 
b. 
= 5.6461 4.5939 4 
5 
8, -.0948 . 2283 
903 
3.6750 3.6681 4 
95; 
2 2% H3484 -.8374 -2 
39, 
202 5.0235 4.9820 5.2060 
25 
IPS) - 75835  .1993 1.1974 
AA 
Ce -.4274  .1673  .7624 
3 
903 
52, 1.1809 ‚3666 1.9598 
903 
= 324 +2, 0191 -o 6132 -3, 3205 
207 
~ f 283 +. 7030 -1.9668 
ESE 
zw! 
AP *04352  * .01116 - 


" % 
1.7698 1.9004 
3.0618 3.2877 
2.0378 2.1757 

‚7433 1.1592 
4.1526 4.7337 
4.0553 4.2569 

.1539 4.1561 4 
‚5665 .9512 1 
.0181 4.7261 5 
.2763 -4.4955 -8 
5.6770 6.4060 
2.3843 3.9077 
1.4158 2.1855 
3.8001 6.0932 
-6.7254 -11.5795 
-5.0229 -9.1524 
.01885 - .05475 - 


2.9912 


3.6175 


2.4178 


1.6730 


5. 8458 


4.8715 


. 5642 


. 4161 


. 8336 


. 2680 


7.4200 


5.9360 


3.1368 


9.0728 


-18, 9730 


-15.2762 


. 10553 
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The constant C, is found by applying Equation (77) at the inner 


channel boundary: 


C, » 3. 7894-. 14286 |a. 8312+2.9234+.0658+2(1. 5221+1.6705)] = 1.7600 
Pz =.5417(.5534)77 =.0683 E (.0683) 7 = .4645 
z 23 b3 c3 ds e3 f3 

5 . 4645 .5080 ^ .5192 .5003 | .4456 . 3400 

A . 0683 . 0934 . 1908 . 0886 . 0591 ‚0229 

Ps 

E ‚1261 ‚1594 ‚1620 ‚1359 ‚0874 .0329 

Tz)" .2278  .2694  .2725  .2404  .1754  .0873 
Average .2486  .2709 .2564  .2079  .1313 

7-1 

1-87) .5355 . 4920 . 4809 . 4996 ‚5543 . 6601 

Average ‚5197 . 4864 . 4902 . 5269 . 6072 


The axial velocity integral is iterated in steps across the channel as 


follows: ad, Pa Average 
,02117 [.3557(5. 1414-=.7702)+3.46=.4073] = .09754 .09754 a 
02125 [.4176(4. 8681-. 8526)+3. 4641. 1633] = .13388 . 23142 ns 
. 02133 [. 4208(4. 9062-. 9333)+3. 46+2. 8802] = .17090 . 40232 ~ 
. 02138 |. 3598(5. 2735-1. 0018)+3. 46+4. 9470|= .21260 . 61492 Im 


.02144 [.2376(6.0766-1.0842)*3.46*7. 5834| - .26220 .87712 

The mass flow — is iterated to solve for Ci as follows: 

. 0292422. 003229(. 0487741. 1799C,)*. 003676(. 16448*1. 0912C,) 
-. 003818(. 31687 41. 0166C, )*. 003522(. 50862 *. 9529C,) 


*. 002257(. 74602 *. 8978C,) 
2.023790 
v, ^TorTOES ^ |. 3920 





Solution at station 3: 


b 


a3 

ph, 1.7098 1. 
2s 2.1980 2. 
O, -.2565 

P, . 0683 

A 1.7601 2. 
Q3 . 03880 

e. 2.8773 3. 
Ost . 13256 

E. . 3815 

M 1.780 l. 
M, 1.100 

AU. K 0001 +, 


10 


3 3 
6727 1.6943 
0458 1.9951 
. 0669 . 3897 
. 0934 ‚1008 
0079 2.1365 
.04653  .04718 
0093 3.1204 
‚12794 .12162 
. 3849 .3795 
570 1.540 
.993 1.005 
0028 -.0170 


d 


3 


1.7698 


. 0378 
. 7433 
. 0886 
. 1406 
. 04139 
3. 
.11358 
. 3641 
. 620 

.110 

. 0548 


2056 


ta 
1.9004 
2.1757 
1.1592 
. 0591 
1.9793 
. 02986 
3. 2572 
. 10373 
. 3379 
1.880 
1. 340 
-. 1065 


3 


. 0912 


2.4178 
1. 6730 


. 0229 
. 5653 
. 01464 
3. 
. 09198 
. 3000 
.420 


2649 


1.810 


. 1603 


Correction of supersonic solution for displacement of streamlines: 


jz 
m/s) "POE d£ 


az B 3 2 e2 = 
$, . 705 . 842 . 899 .956 1.013 1.070 
m(é,) 0 . 005498 .011201 .017071 .023095 .029242 
5 
m(£,} ER Q5 P; 5, dz 
2 
as b3 C3 ds 23 f3 
£, .160 . 848 .936 1.024 lil 19208 
m(g,) © .005258 .011609 .018374 .024639 .029242 


— The functions m( £) and m( $ ) are plotted against radius on 
Fig. 6. From this graph, a corrected ce is obtained and, hence, 


new conditions at station 2. 





Corrected values at station 2: 
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a2 b3 a dz e t$ 
6, 785 . 839 .902 968 1.026 1.070 
Pe .7051 .6778 .6322 .6305 .6157 . 6072 
Ae 1.9390 1.8640 1.7926 1.7338 1.6933 1.6698 
82 -. 0705 . 0090 . 0955 .1804 . 2514 . 3036 
Los 3.7894 3.9438 4.1240 4.3127 4.4785 4.6043 
Py . 5417 . 5837 . 6241 . 6580 . 6816 . 6962 
ie 3.1805 3.3815 3.6043 3.8266 4.0135 4.1519 
O, .17032  .17259  .17315  .17197 .16978  .16767 
Corrected solution at station 3: 

âp b3 c3 a3 2; f3 
P3 1.7398 1.7007 1.7151 1.7686 1.8722 2.0493 
Az 2.2499 2.0942 2.0311 2.0357 2.1269 2.3453 
85 -. 2610 . 1020 . 2573 . 4422 .9361 1.6394 
A . 0605 . 0842 . 0953 . 0959 . 0733 . 0295 
E. 1.7003 1.9450 2.1067 2.2072 2.1223 1.6832 
O5 .03558 | .04330  .04523  .04344 .03452 .01754 
F, 2.8660 2.9896 3.1189 3.2376 3.3064 3.2898 
Dar ‚13104 .12665  .12076  .11325 .10449 .09375 
P4, . 3643 . 3790 . 3763 . 3666 . 3458 . 3071 
AT; +,0006  +.0033 -.0069 -.0368 -.0881 -.1032 
C, = 1.4165 C, = .4487 


1 


2 


Note: Examples VI, VII and VIII use the channel configuration of Fig. 


5 and the same conditions at station l. Example VI is the isen- 
tropic solution (y = 1.40), Example VII is a polytropic solution 
for n = 1.37, and Example VIII is for n = 1.33. From these 
examples it is possible to see clearly the effect of a decrease 
in polytropic efficiency on the solution for the same mass flow, 


upstream conditions and channel configuration. 





A 


e 
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Example Via. 


Solution for the stator of Fig. 5 (divergent curved channel 
walls) with a radial total temperature gradient. This is an isentropic 
(y = 1.40) subsonic solution for the direct problem for a stator 
(F(£,) given). The radial coordinates have been chosen for equal 
areas in the plane perpendicular to the axis. The same equations 


as used in Example Ila are applicable. 


Conditions at station 1: (given) 


ay b, cy di ei fi 
£, .8000 .8438 . 8854 .9252 .9633 1.0000 
T, 3.7894 3.9679 4.1374 4.2995 4.4548 4.6043 
Pig 1.000 1.000 1.000 1.000 1.000 1.000 
Qu .26389 .25202  .24170  .23258 .22448  .21719 
P, .9500 .9500 .9500 .9500 .9500 .9500 
Q, .25439 .24295  .23300 .22421 .21640 .20937 
T 3.7343 3.9102 4.0772 4.2369 4.3900 4.5373 
$, .6202  .6355 . 6481 .6583 .6657 . 6716 
O, -.0310 .0030 . 0367 .0701  .1026 . 1343 
F22.75 m=.026865 2 = 4.0746 A, = 0 
/ 
Solution at station 2: (y = 1.40) 
a, b, C5 d, e f- 
P: .7850 .8496 . 9098 .9661 1.0194 1.0700 
Pe . 5772 .$513 . 5317 . 5166 .5046 .4948 
Ae 1.5873 1.5161 1.4622 1.4206 1.3876 1.3607 
Bo -.897 .0198 . 0865 . 1453 .1986  .2474 
P. .6712 .7095 . 7378 . 7592 .7763 =. 7898 
N, „19848 .19723 .19451  .19103 .18732 .18349 
7, 3.3814 3.5974 3.7930 3.9738 4.1437 4.3041 
M3 .775 .754 .675 . 668 . 630 . 609 
a? 0 -.0013 . -.0025  -.0037  -.0061 -.0080 
e "E. Barton 2 C, = .4661 C, = .8923 





Example VIb 

Solution for the rotor of Fig. 5 (divergent curved channel walls) 
with a radial total temperature gradient. This is an isentropic (y=1.40) 
subsonic solution for the direct problem for the rotor (G( 5.) given). 
The radial coordinates are for equal areas. The same equations as 


used in Example IIIb are applicable. 


Conditions at station 2 are the same as after the stator, Example 


Via. 
G = 1.96 m 2.026865 Fe = 2.8593 
2 
Solution at station 3: (y = 1.40) 
24 "i E Pa e3 fa 
2. . 7600 .8661 .9605 1.0464 1.1258 1.2000 
D; . 1652 . 2370 . 2982 . 3543 . 4089 . 4658 
Ag  -.4362 -.4016  -.3760 -.3520 +-.3244  -.2870 
P4 . 7550 . 7595 . 7626 . 7635 . 7606 . 7532 
Ty 3.4970 3.6683 3.8290 3.9802 4.1197 4.2460 
Oy . 21589 .20707 .19916 .19181 .18460. .17738 
T, 3.5281 3.6993 3.8605 4.0126 4.1549 4.2867 
Pay  -7788 . 7824 . 7847 . 7853 . 7837 . 7788 
Qa, — 22073 21150 .20326 ‚19572 .18859 .18167 
algz  -.0001 -.0002 ~.0024 -.0072 . -.0135  -.0219 
C, = -1329 C, = .9228 


Example Vila 


Solution for the same channel configuration, mass flow, and 
conditions at station l as for Example Vla, except for a polytropic 
process (n 21.37, 727294. 59/0) through the stator. 

E" Equations (48), (57), (58), (66),(74), (75) and the following 


integral are applicable for the polytropic stator (F = const): 











227. «e 
eee E b C 
2 


(Ira) Do AEG, ds, £ LB 
e 
dlie 
= 2.75 = . 026865 =<" = 4,0746 

as 

Solution at station 2: (n = 1. 37) 
a, b, Co d, e, f, 

a .7850 .8496 .9098 .9661 1.0194 1.0700 
Po . 5897 ‚5652 .5469  .5330 ‚5222 ‚5135 
Aa 1.6216 .5543 1.5040 1.4657 1.4360 1.4121 
0, -.0590 . 0203 ‚0890 .1499 . 2055 , 2567 
P, . 6450 . 6837 .7127  .7347 ‚1523 . 7662 
ob 3.3635 3.5778 3.7726 3.9526 4.1218 4.2814 
Q, ..19176 .19110  .18889 .18585 .18251 . 17895 
P. .9789 .9821 .9845  .9863 .9874 . 9882 
Q5, .25834 .24751 .23793 .22938  .22164  .21461 
P *. 0001 -.0017 . -.0022 -.0038  -.0066 -.0090 
Example VIIb 


Solution for the same channel configuration, blade shape, 
and mass flow as for Example VIb, except for a polytropic process 
through the rotor (n = 1. 37, 7 “a 94. 59/o). 

Conditions at station 2 are the same as after the stator, 
Example Vlla. 

Equations (48), (57), (58), (70), (75), (77) and the follow- 
ing integral are applicable for the polytropic rotor (G = const, 


Ps? 1): 





f; 
E / B” £, 27 die Y nd /P 
$ = —_ =— a JEL (È) d ET KAlres 
F 8%, t) | 27 ” = == Lfe à) 
£. 


4 G-a) 9 dar ¿Es 
45 d (72 1)(2-1MX Qe de 


G=1.96 m=.026865 zu z 2.8593 
Solution at station 3: (n = 1. 37) n 

^ a, 95 c3 d3 E f3 
5, . 7600 . 8661 .9605 1.0464 1.1258 1.2000 
P . 1664 . 2496 . 3120 . 3637 . 4097 ,4514 
Az -.4339 . -.3769 +.3490 -.3335 . -.3228 -.3153 
8; -.0250 . 0197 . 0883 , 1704 . 2621 . 3611 
3 . 7387 . 7430 . 7457 . 7463 . 7439 . 7379 
7; 3.4890 3.6590 3.8190 3.9694 4.1094 4.2382 
Q4 .21172  .20306  .19524 .18799  .18101  .17410 
T. 3.5199 3.6839 3.8422 3.9947 4.1404 4.2807 
Ps, . 7620 . 7608 . 7617 . 7632 . 7636 . 1640 
Ost .21646 .20653  .19823 .19102  .18443  .17849 
al; 0 -. 0044 -.0092 -.0137  -.0177  -.0194 
C, = „1338 C, = .9215 


It is interesting to note that, compared to the isentropic case 
(Examples VIa,b), the polytropic flow through the identical stator 


and rotor results in the following percentage changes at station 3 


maximum): 
$4: 55. 39/o | ly -0. 39/o 
Agr LAPIS Py, 3 -3.0%0 


Pa: LA. 3€ 1o Work out: +6. 6?/o 
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Example Villa 


Solution for the same channel configuration, mass flow, and 


conditions at station l as for Example Via, except for a polytropic 


process (n = 1.33, 7 = 86. 8°/o) through the stator. The equations 


used in Example Vila are applicable. 


F = 2.75 m s .026865 a = 4.0746 
/ 
Solution at station 2: (n = 1.33) 
22 Ba .z dg "; f3 
6: .7850 .8496  .9098  .9661 1.0194 1.0700 
Do .6225 .5994  .5828 . 5704 .5613 .5545 
Ào 1.7119 1.6483 1.6027 1.5686 1.5436 1.5249 
0, -.0622 .0216 .0948 . 1605 .2209 .2772 
b. . 5878 . 6280 .6584  .6817 .7000 .7146 
e 3.3149 3.5286 3.7225 3.9022 4.0695 4.2279 
O, .17731 .17796 .17685  .17470 .17199 . 16902 
E». 9387 9470 .9532 9571 .9608 — .9631 
QA .24770 .23863 .23035 .22262 .21564 .20918 
al; -.0001 -.0002 -.0019 -.0044 -.0071  -.0107 
Z Tu CQ 5.5021 C, 2.8765 


Example VIIIb 


Solution for the same channel configuration, blade shape, and 
mass flow as for Example VIb, except fof à polytropic process 
through the rotor (n = 1. 33, 7 = 86. 89/o). 

Conditions at station 2 are the same as after the stator, Ex- 
ample Villa. The equations used in Example VIIb are applicable. 

One compiete approximation for the polytropic, subsonic, direct prob- 


lem for the rotor is given below to demonstrate the iteration procedure. 





Data used in the solution: 





71 


e, 


2 
IE 026065 cedes Zen Zu 
14G 
» dder 
EN aser 210.0075 
74 df, 
a, C d e. f. 
53 .7600 .8661 9505 1.0464 1.1258 1.2000 
af 
E -.1500 .0791 ‚2829 ‚4684 6398 ‚8000 
UE 2500 1720 5474 8891. 1.2049 1.5000 
229, æ o Ca ° e = e . e 
Factor ‚02083 ‚01874 .01716 .01593 .01495 
Te 1.6278 1.7561 2.855) 1.9330 1.9949 
e 72 
S (2 x,-fA) 9.1313 9.0490 8.9897 8.9174 8.8593 
£ 2 
P 
sie | 9.6600 9.8311 10.0646 10.2042 10.2714 
7-1 ds, Qə 
em (B* 2.0529 2.1450 2.2328 2.3169 2.3975 
(%-1)(n-1) \ Qe a 2 . > > 
oe 1.2433 1.1208 1.0325 ‚9646 ‚9102 ‚8653 
5, 
mb e E) 
QUE Es ‚01532 .01530 — .01516 .01495  .01470 
"2 39 ‚006801 ‚006099  .005064  .004509  .004072 


(For definition of "Factor" above, see Problem Vb). 
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ary C, £ .r568 





+; b3 C3 d; d. t 
P; ‚1949 ‚2717 330% „3775 ‚4164 ‚4524 
Go; . 3820 5325 .6492 .7399 .8161 .8867 
Az -.3780  -.3336  -.3113  -.3065  -.3097 -.3133 
8; -.0292 0215 .0937 .1768 .2664 .3619 
n: 
A3 ‚1880 ‚1285 ‚1009 ‚0898 ‚0852 .0818 
És 
Az „1582 ‚1147 ‚0953 .0875 .0835 
5, 
8, -.0039 .0576 ‚1352 .2216 ‚3141 
a pees +.0019 -.0418  -.1308  -.2501  -.4043 
327 
p: 4.3207 -.1296 -.3559 -.4518 -.4635 -.4084 
DN ez), 
(24) -.0487 .0467 ‚1813 ‚3356 ‚5017 ‚6786 
310827 | 
_965 +.2720 -.0829 a re  +.2702 
af, 
265 1.0946 -.1287  -.1454  .-.0390 4.1542 
24, 
Ge 221 4.0388 +.0515 +.0155  -~.0670 
2942 
22-1 


AD x 4.00094 — 4.00068 4.00008  -.00066  -.00158 


The constant C5 is found by applying Equation (77) at the inner 


channel boundary. 


Za = 3.7894». 14286 [ 1429 +. 0380+. 0009+2(1. 3438-. 2873) = 3.4616 


n n-1 
P, = .5878(1. 0443)" = .7000 C, =(.7000) ™ =.9153 
a3 b3 €3 d3 e, f3 
nel 
P,n  .9153  .9162  .9169  .9170  .9163  .9148 
P, ‚2000 .7028 .7049 .7052 .7031 .6984 


P 
(52) 1.1909 1.1191 1.0706 1.0345 1.0044 .9773 
2 


l 

ms 

(52)" 1.1404 1.0883 1.0526 1.0258 1.0033 .9829 
2 


Av. 1.1143 1.0704 1.0392 1.0145 .9931 
n-1 

P, $ .9157 .9165 .9169 .9166 .9155 
d pr 

Lp” +,0380 +.0227 +.0026 -.0209  -.0490 

$, df 


The axial velocity integral is iterated in steps across the channel as 


follows: T s me 
, 02083[1. 8139(9. 1313-8. 8461)+3. 92+. 0946+. 0780] = . 09602 . 09602 : 04901 
. 01874[1. 8797(9. 0490-9. 0107)*3. 92-. 1287 *. 0487] 2 . 07331 .16933 " a 
.01716[1.9278(8.9897-9.2287)+3.92-.1454+.0058] = .05697 . 22630 m 
.01593[1.9610(8. 9174-9. 3537) *3.92-.0390-. 0484] - . 04742 . 27372 "cui 


. 01495[1.9811(8. 8593-9.4040) 3. 924. 1542-. 117512 . 04302 . 31674 


The mass flow equation is iterated to solve for C, as follows: 
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. 026865=. 017071(. 04801 +1. 1820C,)+. 016377(. 13267 +1. 0766C,) 
+, 015754(. 19781+. 9985C,)+. 015167(. 25001 +. 9374C,) 
+, 014597(. 29523+. 8877C,) 


C, - 922474 =. 1568 
Solution at station 3: (n = 1. 33) 
i Ba C3 dy es f3 

E .7600 .8661 .9605 1.0464 1.1258 1.2000 
Ps . 1949 ‚2717 . 3312 ‚3775 ‚4164 „4524 
À3 -.3780  -.3336 -.3113 .3065  -=.3097  -.3133 
0, -.0292 . 0215 . 0937 . 1768 . 2664 .3619 
P3 . 7000 . 7028 . 7049 . 7052 . 7031 . 6984 
De 3.4617 3.6285 3.7862 3.9354 4.0740 4.2038 
Qz .20220  .19367 .18615 . .17921  .17275  .16518 
la. 3.4876 3.6503 3.8062 3.9582 4.1049 4.2456 
P4. . 7185 7147 . 7181 . 7204 . 7219 . 7229 
Qs, .20600 .19659  .18863 .18183  .17605  .16929 
AT; -. 0001 -.0047 -.0108 -.0155  -.0178 -. 0202 


Compared to the isentropic case (Examples Vla,b), the poly- 
tropic flow with n = 1,33 through the same stator and rotor results in 


the following maximum percentage changes at station 3: 


$4: +18. 0°/o un. -1. 19/o 
Ay: +20. 8°/o Pat -9. 39/o 
P4: - 8. 39/o Work out; 119. 59/o 


Example Lia 


Solution for the stator of Fig.5 (divergent curved channel 
walls) with a radial total temperature gradient and radial coordinates 
for equal areas. This is an isentropic, subsonic solution for the direct 


problem for the stator. The resulting solution at station 2 is used in 





Bi 


zxample [Xb as the upstream conditions for a solution of the inverse 
problem for a rotor. 


The equations used in Example Ila are applicable. 
Conditions at station 1: given 


a 5 e di e A 
5, .8000  .8438  .8854  .9252  .9633 1.0000 
B 3.7894 3.9679 4.1374 4.2995 4.4548 4.6043 
P. 1.00 1.00 1.00 1.00 1. 00 1.00 
Qu .26389  .25202  .24170  .23258  .22448  .21719 
P, .9400 . 9400 . 9400 . 9400 .9400 . 9400 
i, 3.7231 3.8985 4.0650 4.2243 4.3768 4.5237 
Q .25249  .24113 .23126 .22253 .21478 .20781 
8, -. 0340 . 0033 . 0403 . 0768 . 1125 .1473 
$, . 6804 . 6970 . 7108 .7215 . 7303 . 7365 
F 22.75 m 2.029242 ig z 4. 0746 À=0 

/ 

Solution at station 2: 

a: b2 e? d, ez E 
A . 7850 . 8496 . 9098 .9661 1.0194 1.0700 
Po . 7008 . 6688 . 6452 . 6275 . 6139 . 6033 
Az 1.9272 1.8392 1.7743 1.7256 1.6882 1.6591 
0; -.0701 . 0241 . 1050 . 1765 . 2416 . 3016 
P, . 5460 .5954 . 6322 . 6601 . 6822 . 6997 
=, 3.1879 3.4217 3.6292 3.8183 3.9938 4.1577 
O, .17126  .17402  .17419  .17286  .17081  .16828 
alo. 0 -.0010  -.0026  -.0049  -,0083  -.0116 

PET P»? Pug Soe = Me 
C, 2.5659 C4 2.8412 

Example IXb 


Solution for the rotor of Fig. 5 (divergent curved channel 


walls) with a radial temperature gradient and radial coordinates for 





equal areas. This is an isentropic, subsonic solution for the in- 
verse problem (H(£. ) given) for the rotor. One complete approxi- 
mation for the inverse problem for the rotor is given below to il- 
lustrate the iteration procedure. | 
Equations (63), (65), (73), (75), (77) (with appropriate 


subscripts) and the following two relations are used: 





Psy £s | v atio, E ees) dca Au 
E 2/27 S dE) * dg dle, a ie dE aeo (80) 


H(£)7 £ À,* £ A, (27). 


When iterated across the channel, Equation (80) becomes: 





2.22 E tes js PEE fi) (e t A A gat , 
(81). 


* a (H 4.) v Cn te) 


Conditions at station 2 are the same as after the stator, 
Example IXa above. The work output function E(4, ) is given in thie 
example, and the blade shape (spouting angle) is determined. 


Data used in the iteration: 





Aloe . » (constant- 
= ,029242 qz, "2.8593 H(f,) 2 1.300 ver) 
dH 


(Mi-am 0 a£ 
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23 b3 C3 d, e f3 
Z .7600 .8661 .9605 1.0464 1.1258 1.2000 
H(é,) 1.3000 1.3000 1.3000 1.3000 1.3000 1.3000 
TUR 1.5129 1.5626 1.6143 1.6671 1.7210 1.7752 
mA, -.2129 -:32826  -.3148 -.3671 24210 -.4752 
Az -. 8801 +.3032 -.8272 -.3508 -.3739 =. 3960 
mU 
E 1.7889 1.7317 1.6766 1.6238 1.5736 
nes 
2 -e 
z di £i) .86305 .86305  .86305  .86305 .86305 
ds, e 
(H- PA, berry fe -.02931 -.02801  -.02751 -.02681 -. 02631 
2 p 
E.) .1061  .0944 .0859 .0794 .0742 
= (£- £5) 
a ‚01489 .01501  .01496 .01482 .01462 
.007820 .006658 .005821 .005185 . 004682 


Z Ef) 





Try Ci =.1039. 
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ay b3 c3 “ds e3 t 
P; , 1039 ,2379  .3295  .4159 ‚5242 .6891 
Az .. 2801  -.3082  -.3272 ~.3508  -.3739  -.3960 
03 -. 0156 . 0188 . 0932 ‚1948 ‚3354 ‚5513 
I 
A3 
. 1047 , 1088 ‚1145 ‚1209 ‚1274 
53 
O, ‚0016 ‚0560 . 1440 .2651 ‚4433 
aS ‚3242 ‚7881 1.1827 1.7707 2.9097 
5, 
- 3,22 -.0005  -.0441  -.1703  -.4694 -1.2889 
2 3 
36 -2,9845 -2.1545 -1.5785 -1.0580 -.4485 +.4290 
DPIBE) 4477  -.1704 -. 4466  -.4956  -.2870  +.3482 
2 3 ec 7; 
" | 
$575) -.0260 ‚0409 .1804:  .3698 ‚6316 1.0336 
Jig 7 
203 +.4217  -.1295 -.2662  -.1258  +.3446 +1.3768 
262 
903 +1461 — -.1978  -.1960 3.1094  +,8607 
2c, 
- ¿20% -.0065  +.0592  +,0700 -.0642  -.5647 
226 
Z 
AE +,00076 +.00082  :.00008  -.00214 -.00809 
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* 
The constant C, is found by applying Equation (77) at the inner 
channel boundary. 


T; = 3.7894 - .14286 [.0785 + .0108 +. 0002 + 2.600] = 3.4052 


2 7 
Pz =.5460 (1.0682)7=7 -.6878 GC, =(.6878)7 =.8986 
a3 b3 cz d3 ez fy 
7-1 
Pj"  .8986 . 8994 ‚9002 ‚9003 ‚8981 ‚2900 
P, . 6878 . 6900 ‚6921 . 6924 ‚6865 , 66 
(E) 1.2597 1.1589 1.0947 1.0489 1.0063 9506 
zb 
(-2)" 1.1793 1.1111 1.0668 1.0347 1.0045 . 9645 
e 
— 1.1452 1.0839 1.0507 1.0196 .9845 
TA 
(1-8 7).1014 ‚1006 . 0998 .0997 ‚1019 ‚1100 
E. ‚1010 ‚1002 0998 . 1008 . 1059 


The axial velocity integral is iterated in steps across the 


channel as follows: 


AP Py Average 
2.0486 (.08717 .02931) + .01550 = .13403 .ı3403 06701 
1.8856 (. 08648 -. 02801) -.01867 =.09158 . 22561 17782 
1.7616 (. 08613 -. 02751) -. 01684 = .08642 .aı203 26992 
1.6556 (.08700 -.02681) *.00869 2.10834 .42037 p 


1.5492 (.09140 - .02631) +.06386 = .10470 . 58507 


The mass flow equation is iterated to solve for Ci ag foilows: 


c 


2)242 2 .017052 (. 06701 * C,) * . 016344 (. 17982 * C) 
+ . 015718 (. 26882 + C,) * .015110 (. 36620 * C) 
+. 014393 (. 50272 + C,) 


008166 


5 = .1039 


GQ - 
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Solution at station 3: 


a3 ag cg d3 

$,  .7600 .8661 .9605 1.0464 
$,  .1039 . 2319 . 3295 ‚4159 
Az -.2801 -.3032 -.3272 -.3508 
0, -. 0156 . 0188 . 0932 . 1948 
P .6878 . 6900 .6921 . 6924 
T, 3.4052 3.5688 3.7243 3.8766 
O4  .20197  .19335  .18583 .17886 
Gf) 4.619 2. 366 1.922 1.672 
B, 12 12? 22558:  27%30'  30%55' 
aT; 0 1.0064  +.0096  +.0097 


“E 1; 
1.1258 1.2000 
. 5243 . 6890 
-.3739  -.3960 
. 3354 ‚5512 
. 6865 . 6651 
4.0010 4. 0978 
‚17158 . 16231 
1. 434 1,167 
34°53' 40°35! 
+,0070 — 4.0014 


Note: A supersonic solution is indicated for problem IXb by using 


average values in the mass flow and energy equations 


( $ = 3.160). 


iowever, the supersonic solution is im- 


possible in this case since the pressure goes to zero at the 


outer channel boundary before the solution is reached. 


Exam ple Xa. 


Examples Xa, b are the isentropic, subsonic solutions for the 


stator and rotor of Fig. 5. 


in which H( e ) and 


both the stator and rotor are thus determined. 


This is an example of the inverse problem 


A = 0 are specified and the blade shapes of 


Equations (27), (63), (65), (73), (74), (75) (with appropriate 


subscripts) are applicable for the stator in addition to the following 


relation: 


(A; = O, H($,) - £ À,) 











2 
aA 7 £, xy al OS As Al: 2Az2 
4 - As zs (E IARE nto 0n 


Conditions at station l are assumed the same as for Example 


IXa. 

az by Ez MB: E £ 
H(%) 1.3000 1.3561 1.4061 1.4515 1.4935 1.5328 
(given) 


1.6561 1.5962 1.5455 1. 5024 1.4651 1.4325 





m = .029242 oe 
gon = 4. 0746 Fir Pe 
7 
Solution at station 2: 
à b, C5 d, e, f, 

6,  .7850 .8496 .9098  .9661 1.0194 1.0700 
pə .6791 . 6327 . 6004 . 5760 . 5532 . 5311 
Afevw)). 6561 1.5962 1.5455 1.5024 1.4651 1.4325 
O,  -. 0679 . 0228 . 0977 ‚1624 ‚2177 . 2656 
D, «6867 . 6180 . 6932 . 7308 '. 7489 . 1635 
T, 3.3310 3.550] 3.485 3.9 mic 4.1015 4.2627 
Q,  .19116  .19093  .18604  .18590 . 18279  .17912 
F(€,) 2.439 2. 52$ 2. 574 2. 608 2. 648 2. 697 
Be 2218:  21%36:  21%15:  20%59:  20%1'  20%21' 
AT, 0 -.0035 -.0052 -.0051 -.0038 -.0019 


C,=.6791 Cie Ae e We 


e 2t lt 12, 
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Example Xb. 

Solution for the rotor of Fig. 5 for the inverse problem in 
which H ( 5, j and 2 O are specified and the blade shape G (€) 
is to be determined. 

Equations (27), (63), (65), (73), (75), (77) (with appropriate 


subscripts are applicable for the rotor in addition to the following 


p 
214 ef of! 
een © 
$; 


Conditions at station 2 are the same as after the stator, 


Example Xa. 


= , 029242 SF = 2.8593 A,= 0 Aé, ef 

Solution at station 3: 
a3 ES cs d3 tg is 

£, «7600 .8661 .9605 1.0464 1.1258 1.2900 
ds .2985 . 3398 . 3627 . 3830 .4110 . 4602 
83 -.0048 . 0269 . 1026 .1794 . 2630 . 3682 
Az! 0 0 0 0 0 0 
P,  .6876 . 6873 . 6868 . 6857 . 6825 . 6758 


Q, 20196 .19280 .18481 .17763 .17105 . 16416 
T, 3.4049 3.5639 3.7388 3.8602 3.9940 4.1165 
G($,) 2.546 2.549 2.648 2.732 2.739 2.608 

Q; 21°25" 210723 Bm aor 20°02" 207 MP. 
al, 0 -.0001 4.0011 -.0010 +.0001 +.0002 


C, =.2985 C, = .8985 


l 2 
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Example XI 

This example has been constructed to illustrate the increase 
in gross thrust obtainable by the use of a radial total temperature 
gradient. 

A hypothetical one-stace turbine consisting of a stator and a 
rotor has been assumed similar to that shown in Fig. 3. In order 
to give this example more physical meaning, rotor blades have been 
used which are tapered linearly in area in such a way that the ratio 
of hub to tip cross-sectional area is three. Using a rotor tip speed 
of 1100 feet per second, the stress due to centrifugal forces was 
calculated at various pointe along the blade. Using the high-tempera- 
ture data for 'Vitallium'', (H. 5. 21) (NR-10), given in References 16 
and 17, it was then estimated that the blade root could safely stand a 
temperature of 1450°F, while the blade tip temperature could be 
1850°F. 

Using this data, the thrust of this turbine is calculated for 
two conditions: (a) A constant radial total temperature of 1450°F í 
and {b) a linear temperature gradient of 400°F with the same root 
total temperature of 1450°F. The mass flow, boundary configuration, 
combustion chamber total pressure, and wor extracted arc the same 


in both cases. 


Example „la. 


Solution for the stator and rotor of Fig. 3 with no radial total 


‘temperature gradient and equal area coordinates. This is an isen- 


tropic, subsonic solution to the inverse problem (H /( ay, given). 
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. 
For no radial total temperature gradient, H ( = )= A 
= constant, and a cylindrical channel, Equation (82) reduces to 
BP, = C,- Equations (63), (64), (65) and (67) are also used in 
the solution for the stator. 


Conditions at station 1: (given) 








222.7071 Pj51.00 — QC, - .36940 P, = .4915 
) =26726  — PF, - .9962 Q, = .35775 A,= 0 
3 rad E O m lb 
w = 880 rie TE 1910"R (const.) W = 148.26 E 
wR= 1100 E pe LO, O H(%)= 1.3964 (const .) 
o sec P ! «^ 3 " ; 
d Oft m = .031650 Àz=0 
Solution at station 2: 
a, b, c d, | e, 
3 .8000 ‚8544 ‚9055 ‚9539 1.0000 
| Po ‚7260 .7260 ‚7260 .7260 ‚7260 
| A> — 1.7455 1.6342 1.5421 1.4639 1.3964 
1 .4812 5236 ‚5585 ‚5879 ‚6125 
e. .21909 .23268 .24369 ‚25277 26027 
t 2.1965 2.2501 2.2921 2.3259 2.3533 
F - 4 2.404 2.251 2.124 2.016 1.923 
2 
Pe 22934: 23°59) 25°13" 26925: 2358 39! 
M, 1.071 1.605 ‚951 .904 .873 
ato 0 +.0002 0 - .0003 -.0001 
C} = -7620 la ee ae ee ee 
C, - .8114 


For the rotor with no radial total temperature gradient, 
— E (£, ) - constant, end A - 0, the solution at station 5 


is easily found to be: 


Pz = .7806 Q4 7 .22525 
A de = 
3 0 Lo. =2- 3081 

D = £ 

P. . 5003 = . 24798 
fm = 2.2211 Pa, = -5724 

az b3 C3 d3 *3 
G = i 1.025 1.095 1.160 1.222 1.281 


85 44?20' 42?25' 40°46' 39°17' 37756 
Assurning isentropic expansion to atmospheric pressure through 


an ideal nozzle, the total gross thrust is found to be 11, 966 pounds. 
Example Älb. 

Solution for the same stator and rotor as Example Xla except 
with a linear radial total temperature gradient of 4007. The mass 
flow, combustion chamber total pressure, and work extracted are the 
same as in the above example. 


Equations (27), (63), (64), (65), (67), and (82) are applicable. 


Conditions at station 1: (given) 


a] E S x £j 

> . 8000 .8544 .9055 .9539 1. 0000 
p. 2. 0871 2.8614 3.0062 3.1434 3.2741 
Pit 1.00 1.00 1.00 1.00 1.00 
Ort . 36940 . 34943 .33265 .31813 . 30543 
Pi .9512 .9512 .9512 .9512 .9512 
T, 2.6687 2.8208 2.9635 3.0988 3. 2276 
a . 35643 .33721 .32097 .30696 , 29471 
$, . 5185 ,5331  .5467 . 5583 .5705 


A, O 0 0 0 0 
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, 
m 7.031650 mA = 2.8350 H(%) 1.3964 (const.), À* C 
Solution at station 2: 
= Ba S5 e *4 
Pe . 5844 .6975 | . 7955 . 8819 . 9589 
A2 1.7455 1.6342 1.5421 1. 4639 1. 3964 
P, 5018 .5443 .5774 . 6036 . 6249 
© 20831 24050 2.5697 27212 2.8625 
C 22574 .22632 .22474 .22183 . 21829 
E de 2. 987 2.342 1.938 1.659 1.467 
Bo 18°30! 230004 27918: año 34719! 
M, 1.044 .974 .920 . 878 . 850 
a1; -.0001 +.0054 +. 0064 +0049  +.0017 
Sr aa ee". + Pa" Pre Sot = Oe 
C, = -8212 


Tor the rotor with E (=) = constant and A; = 0, 
Equation (64) reduces to Pg 7 constant. Equations (63), (65), (71), 


and the following integral are applicable: 





i, 
ANN. 7 dl = 
% 2 qe vu (|-8* Mf c, (84), 





Solution at station 3: 


3 3 

P3 6002 .7426 

Az O 0 

F, 5289 5289 

T, 2.2566 2.3853 

Q, 23435 .22171 

m. 2.3081 2.4624 

P. 5724 ‚592 

Fr 24798 .24005 

Boe). 333 1.151 

Ps 

E 36^ 52: 40958' 

az, 0 8617 
C, - .6002 


.6514 
0 


.5289 


2.5060 


21103 


2.6072 


.6075 


.23300 


1.064 
43°13! 


- .0023 


9397 
J 


‚5289 


2.6203 


.20182 


2.7444 


‚22658 


1.015 
44°35! 


- .0020 


= .8336 


1.0141 
0 

‚5289 

2.7293 
.19376 

2.8751 
6345 
.22067 
.986 
45°22! 


4.0011] 


Assuming isentropic expansion to atmospheric pressure 


through an ideai nozzle, and summing the thrust by stream tubes, 


the total gross thrust is found to be 12,883.68 pounds. 


This is 


an increase of 917.8 pounds, or 7.7 percent, over than found in 


Example Xla. 
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FIG.3 — CHANNEL CONFIGURATION AND COORDINATES FOR 
EXAMPLES I, IL, AND XI 
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O 0.10 0.20 0.50 
5 


Axial Goordinate 


FIG.5 — CHANNEL CONFIGURATION AND COORDINATES FOR 
EXAMPLES IY - .X- 





Examples IV, V 


le | eso | .009 | .9s6 
“7 | 
: 





.800 
.840 
.880 
| 
1.070 


Examples VI - X 


ix) 89) (E) (E) (E) 


.0360|'.0791| .1332| .1467| .1720. 
C| .8854| .9c98| .9605| .o567| .1627| .2829. .4972| .5163| .5474. 
á| .o252| .9561|1.0464 | .1cós| .2813| .4684| .8455| .8620| .8891 
e| .96331.0194]1.1258| .1541| .3935. .6398|1.1789]| 1.1893|1.7049 
1.0700 









TABLE T — STREAMSURFAGE COORDINATES, SLOPES AND CURVATURES 
USED FOR EXAMPLES IY. --X. AND FIGURE 5 





Mass Flow 






Station 2,Prob.IVo 
(Subsonic) 







tation3,Prob. V. b 
(Supersonic) 












HU Station 3, Prob.IV b 
(Subsonic) 
wan 


EHE 
d q 


0.76 0.84 0.92 IOO 
Radial LE 


FIG.6 — MASS FLOW FUNCTION FOR SUBSONIC AND SUPERSONIC 
ROTOR FROM SUBSONIC STATOR 
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FIG. 7 —STREAMLINE SHIFT THROUGH ROTOR 
FOR SUBSONIC AND SUPERSONIC 
“GOVERNING MAGH NUMBER 
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FIG.8 —TYPICAL TURBINE VELOCITY DIAGRAM SHOWING SIGN 
CONVENTION USED IN ALL EXAMPLES 
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FIG.9 — TYPICAL COMPRESSOR VELOCITY DIAGRAM SHOWING 
SIGN CONVENTION 
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